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Abstract 

A  finite  difference  computer  model  is  developed  to  simulate  the  exposure  statistics  of  a  radio 
frequency  buoyant  antenna  as  it  is  towed  in  a  three-dimensional  random  seaway.  The  model 
allows  the  user  to  prescribe  antenna  properties  (length,  diameter,  density,  etc.),  sea  conditions 
(significant  wave  height,  development  of  sea),  tow  angle,  and  tow  speed.  The  model  then 
simulates  the  antenna-sea  interaction  for  the  desired  duration  to  collect  statistics  relating  to 
antenna  performance.  The  model  provides  design  engineers  with  a  tool  to  predict  antenna 
performance  trends,  and  to  conduct  design  tradeoff  studies.  The  floating  antenna  envisioned 
is  for  use  by  a  submarine  operating  at  modest  speed  and  depth. 
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Chapter  1 


Introduction 

1.1  Motivation 

As  the  United  States  Navy  launches  the  last  of  the  highly  capable  Seawolf  class  submarines  and 
begins  construction  on  the  even  more  capable  Virginia  class  submarines,  it  continues  to  find  itself 
operationally  limited  by  the  technological  barriers  present  with  high  data  rate  communications. 
The  capability  for  high  speed  data  transfer  has  escalated  at  an  astonishing  rate  over  the  last 
few  years;  however,  the  submarine  force  faces  unique  challenges  in  fully  utilizing  high  data 
rate  systems  which  operate  at  higher  frequencies.  High  frequency  (UHF,  SHF,  or  EHF  band) 
radio  signals  are  unable  to  penetrate  below  the  surface  of  the  water;  so,  while  the  submarine 
is  submerged,  radio  communications  systems  are  limited  to  low  data  rate  receive-only  modes 
transmitted  at  very  low  frequencies  (ELF,  or  VLF  band).  The  current  mode  of  operation  calls 
for  the  submarine  to  routinely  come  to  periscope  depth  to  raise  a  radio  mast  above  the  surface  of 
the  water,  send  and  receive  the  necessary  message  traffic,  then  submerge  to  continue  its  mission. 
This  required  visit  to  periscope  depth  not  only  causes  the  submarine  to  suspend  execution 
of  its  mission,  but  also  reduces  the  overall  stealth  of  the  mission.  Submarine  operations 
are  typically  driven  by  the  need  for  stealth,  and  stealth  is  best  preserved  by  remaining  fully 
submerged.  Clearly,  the  communications  goal  of  the  submarine  force  is  to  be  able  to  utilize 
the  high  frequency  and  high  data  rate  communication  circuits  while  retaining  the  stealth  of 
operation  at  speed  and  depth. 
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1.2  Background 


This  goal  of  high  data  rate  communication  at  speed  and  depth  is  the  motivation  of  a  U.S.  Navy 
Advanced  Technology  Demonstration  (ATD).  Initially  funded  jointly  by  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  and  the  Office  of  Naval  Research  (ONR)  in  FY  96-97  as 
a  feasibility  study,  the  Buoyant  Cable  Array  Antenna  (BCAA)  research  became  a  DARPA  and 
ONR  Demonstration  in  FY  98-99.  The  BCAA  prototype  satisfactorily  demonstrated  receive- 
only  capability  with  a  field  test  in  the  Summer  of  1999.  The  program  was  promoted  to  the  level 
of  Navy  ATD  in  FY  00  and  renamed  the  Multi-element  Buoyant  Cable  Antenna  (MBCA).  The 
MBCA  is  a  3  year  joint  program  which  is  a  combined  Navy  ATD-funded  and  DARPA-funded 
effort.  This  ATD,  one  of  only  6  ATD’s  funded  by  the  Navy  for  fiscal  year  2000,  reflects  the 
high  priority  the  Navy  places  on  submarine  connectivity.  The  goals  of  the  MBCA  program 
include: 

•  demonstrate  24  kbps  2-way  communication  at  UHF  over  the  Fleet  Satellite  (FLTSAT) 
communications  system  while  the  submarine  is  at  a  speed  of  6  knots  and  at  a  depth  of 
300  feet 

•  demonstrate  1.5  Mbps  data  rates  over  UHF  Line-of-Site 

•  demonstrate  L-band  voice  communications. 

The  participants  in  the  program  include  the  Naval  Undersea  Warfare  Center,  Johns  Hopkins 
University  Applied  Physics  Laboratory,  Sippican  Inc.,  and  MIT’s  Lincoln  Laboratory.  Lincoln 
Laboratory  is  a  federally  funded  laboratory  administered  by  the  Massachusetts  Institute  of 
Technology.  Lincoln  Laboratory,  whose  specialties  include  radar  and  electrical  engineering, 
was  chosen  to  overcome  the  technological  challenge  of  designing  a  floating  linear  array  radio 
frequency  communication  antenna  module  and  the  associated  inboard  electronics.  However,  to 
gain  insight  into  the  behavior  of  the  antenna  module  in  a  random  seaway,  Lincoln  Laboratory 
added  Professor  Jerome  Milgram  of  the  MIT  Ocean  Engineering  department  to  the  antenna 
development  effort.  Professor  Milgram  joined  a  working  group  that  was  established  to  estimate 
statistics  on 

•  the  average  fraction  of  antenna  elements  exposed  over  an  extended  period  of  time 
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•  the  temporal  variation  of  the  exposure  statistics,  including  duration  of  reception  dropout 
times 

•  the  spatial  distribution  of  those  exposed  elements. 

This  thesis  work  was  undertaken  to  support  the  statistics  working  group  by  improving  upon 
the  first  hydrodynamic  model  which  was  developed  by  Professor  Milgram  and  Gary  Ulrich. 
[9]  [10]  The  improved  model  retains  all  of  the  capabilities  of  the  original  model  and  improves 
upon  the  original  model  by  upgrading  the  wave  field  from  a  unidirectional  wave  field  with  the 
choice  of  only  two  tow  directions,  0  or  180  degrees,  to  a  multidirectional  wave  field  with  the 
choice  of  any  tow  direction  from  0  to  180  degrees.  This  provides  all  angles  due  to  port  and 
starboard  symmetry. 

1.3  Antenna  Model 

The  method  used  to  determine  the  antenna  behavior  and  the  exposure  statistics  was  to  develop 
long  term  simulations  from  which  exposure  statistics  could  be  obtained  vice  getting  the  exposure 
statistics  from  the  theory  of  a  non-linear  process.  The  model,  named  PROTEUS ,  performs 
the  simulation  by  constructing  a  wave  field  and  then  modeling  the  antenna  in  the  wave  field 
and  analyzing  the  resultant  behavior.  The  development  of  the  wave  field  model  is  described  in 
Chapter  2  and  the  development  of  the  antenna  model  is  described  in  Chapter  3.  The  results 
of  the  improved  model  are  presented  in  Chapter  4  and  a  detailed  user’s  guide  is  included  as 
Appendix  A. 
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Chapter  2 


Wave  Field  Model  Development 

2.1  Linear,  Random,  Deep  Water  Wave  Theory 

To  simplify  the  model  and  to  reduce  computation  time,  the  simulated  wave  field  was  assumed  to 
be  described  by  linear,  random,  deep-water  wave  theory.  The  linear  theory  aspect  of  the  waves 
implies  that  the  velocity  potential  and  fluid  velocity  are  constant  from  the  mean  free-surface  to 
the  free-surface  level.  The  linear  theory  assumptions  represent  a  first  order  approximation  in 
satisfying  the  free  surface  kinematic  and  dynamic  boundary  conditions.  A  complete  derivation 
and  discussion  of  linear  wave  theory  can  be  found  in  Faltinsen  and  Newman.  [4]  [6]  “Deep 
water”  is  defined  to  exist  when  A  <  2d,  where  A  is  the  wavelength  and  d  is  the  water  depth. 
The  random  aspect  of  the  waves  implies  that 

•  the  instantaneous  wave  elevation  from  the  mean  free-surface  is  gaussian  distributed  with 
zero  mean  and  variance 

•  the  waves  are  an  ergodic  process,  which  means  that  the  wave  statistics  obtained  from  a 
single  location  are  the  same  as  the  wave  statistics  obtained  from  an  ensemble  of  locations 

•  the  waves  are  a  stationary  process,  which  means  that  the  wave  statistics  are  invariant 
with  respect  to  time  [7]  [8] 

Because  the  last  condition  is  so  restrictive,  we  relax  it  to  allow  the  statistics  to  vary  slowly 
with  time,  otherwise  known  as  a  weakly-stationary  process.  A  summary  of  the  deep  water 
linear  theory  wave  relations  is  provided  in  Table  2.1. 
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Quantity 

Formula 

Wave  Number 

k=j 

Circular  Frequency 

Dispersion  Relation 

oj  =  y1 

k  =  T 

Velocity  Potential 

<f>  =  ^fekz  sin(kx'  —  ut) 

Surface  Displacement 

h  =  acos(kxf  —  u ut) 

Radian  Frequency 

u  =  yfgk 

Phase  Velocity 

c=  y/g/k 

Group  Velocity 

Cg  =  \\fgJk 

Horizontal  Fluid  Velocity 

u  —  au>ekz  cos  (kx1  —  ut) 

Vertical  Fluid  Velocity 

w  =  auekz  sin(fcE,  —  ut) 

Wave  Energy  Density 

E  =  \pga 2 

Note:  xf  refers  to  a  stationary 

x  coordinate  system. 

Table  2.1:  Linear  Theory  Deep  Water  Wave  Relations 


Assuming  a  lower  wave  frequency  bound  of  0.3  rad/ sec,  which  is  reasonable  for  the  vast 
majority  of  energy  density  spectra  one  would  encounter,  this  translates  into  a  minimum  water 
depth  of  342  meters.  This  implies  that  for  depths  shallower  than  about  350  meters,  the 
deep  water  theory  begins  to  break  down,  and  a  finite  depth  treatment  of  the  waves  must  be 
considered.  A  complete  derivation  and  discussion  of  finite  depth  wave  theory  can  be  found  in 
Faltinsen  and  Newman.  [4]  [6] 

Upon  inspection  of  the  linear  theory  formulas  in  Table  2.1,  the  waves  are  found  to  be 
sinusoidal.  However,  an  observation  of  the  sea  surface  shows  that  the  waves  display  a  quite 
complicated  and  random  structure.  This  observed  behavior  can  be  numerically  approximated 
by  adding  together  a  series  of  sinusoids  with  each  discrete  frequency  component  having  the 
appropriate  amplitude  and  a  random  phase.  As  the  number  of  discrete  frequency  components 
increases,  the  numerical  approximation  approaches  the  continuous  distribution  of  the  actual 
sea  surface.  This  discrete  method  of  describing  a  two  dimensional  random  seaway  is  presented 
in  Equation  2.1 


N 

h(x',  t )  =  ^  ai(uj)  cos(ki(u)x'  —  u^t  +  <p/)  (2.1) 

i—1 


where 
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Bretschneider  Energy  Density  Spectrum 


Figure  2-1:  Typical  Wave  Energy  Density  Spectrum 

h(x',  t )  is  the  wave  elevation,  referenced  to  the  mean  free-surface 

ai(oj)  is  the  appropriate  wave  amplitude  of  each  discrete  frequency  component 

ki(u>)  is  the  wave  number  of  each  discrete  frequency  component 

Ui  is  the  set  of  discrete  circular  frequencies  contained  in  the  modeled  energy  density  spectrum 
fa  is  the  random  phase  assigned  to  each  discrete  frequency  component 
and  the  ekz  term  has  been  dropped,  since  evaluation  occurs  at  the  sea  surface  (z  =  0). 
Using  the  equations  found  in  Table  2.1,  one  can  determine  most  of  the  terms  on  the  right 
hand  side  of  Equation  2.1:  k{  is  determined  from  the  dispersion  relation,  and  fa  is  a  random 
number  with  a  uniform  distribution  on  0  to  2-n.  The  determination  of  the  appropriate  wave 
amplitude,  however,  requires  some  discussion  and  a  deeper  understanding  of  wave  theory. 

To  determine  the  wave  amplitude  corresponding  to  a  given  wave  frequency,  one  begins  by 
defining  the  spectral  energy  density  function,  S(oj),  which  represents  the  energy  of  the  waves 
having  frequency  u.  A  typical  energy  density  spectrum  is  shown  in  Figure  2-1.  The  units  of 
S(lj)  are  defined  to  be  in  m2-sec,  so  that  when  S(uj)  is  integrated  with  respect  to  uj,  the  result 
is  the  average  energy  of  the  random  waves  with  respect  to  time. 
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Bretschneider  Energy  Density  Spectrum 


Figure  2-2:  Determination  of  Component  Amplitudes  Oi  (u>) 

The  spectrum  shown  in  Figure  2-1  is  said  to  be  a  single-sided  spectrum,  which  refers  to 
the  fact  that  the  frequencies  are  constrained  to  be  positive.  For  a  double-sided  spectrum,  both 
the  positive  and  the  negative  frequencies  are  allowed,  and  the  spectrum  is  defined  such  that 
S(uj)  =  S(—lu).  In  addition,  the  double-sided  spectrum  is  half  the  magnitude  of  the  single-sided 
spectrum;  since,  for  constant  sea  severity,  the  area  under  the  spectral  energy  density  curve  must 
be  constant.  The  amplitudes,  ai(uj),  of  the  discrete  frequency  components  in  Equation  2.1  are 
obtained  from  Equation  2.2 

iaiH2  =  S(w)Aw  (2-2) 


where 

o*( u)  is  the  amplitude  of  the  discrete  band,  Au,  of  waves  with  center  frequency  u. 
Equation  2.2  is  simply  a  statement  that  the  square  of  the  wave  amplitude  is  a  measure  of 
the  average  energy  of  the  waves,  because  both  sides  of  the  equation  represent  energy  in  units 
of  m2.  The  method  of  obtaining  the  amplitudes  cq( u)  is  shown  in  Figure  2-2. 
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Until  now,  S(u>)  has  been  presented  in  terms  of  wave  energy  present  at  a  given  wave  fre¬ 
quency.  It  can  be  shown  from  the  Wiener-Khintchine  theorem  that  S(uj)  is  related  to  the  wave 
amplitude  record  at  a  fixed  location  in  the  ocean.  [7]  We  define  a  wave  record  at  a  given 
location  as  Then  the  autocorrelation  of  xf(t)  is  given  by 

R(t)  =  lim  ~  f  x'{t)x'(t  +  T)dt  (2.3) 

1  —*00  Zi  J  —T 

and  by  the  Wiener-Khintchine  theorem,  the  spectral  energy  density  function,  S(u),  and  the 
autocorrelation  function,  f?(r),  form  a  Fourier  transform  pair. 


sw)  = 

1 

7T  j 

poo 

/  R(t)  e~™TdT 

' — oo 

(2.4) 

R(  t)  = 

1 

2  J 

POO 

1  S(uj)  e-^du} 

—  OO 

(2.5) 

This  method  proves  to  be  invaluable  as  it  has  the  ability  to  convert  observed  wave  records  into 
spectral  energy  density  functions. 


2.2  Energy  Density  Spectra  Contained  in  the  Model 

Numerous  energy  density  spectral  forms  have  been  proposed  by  researchers  attempting  to  char¬ 
acterize  a  particular  sea  wave  environment.  [4]  [5]  [7]  [8]  The  spectral  forms  may  characterize  a 
particular  geographic  wave  environment  or  they  may  characterize  a  particular  storm  environ¬ 
ment,  but  in  all  cases,  they  share  a  very  similar  form.  The  model  incorporates  the  spectral 
formulation  as  a  subroutine;  therefore,  it  is  completely  arbitrary  from  the  model  viewpoint  how 
the  spectrum  is  defined.  In  order  to  provide  some  variability  in  the  allowable  spectral  forms 
and  to  allow  for  troubleshooting  and  debugging  of  the  model,  four  wave  spectral  energy  forms 
were  included  in  the  model.  Two  of  the  spectral  forms  are  actual  derived  models  that  are  well 
accepted  in  the  ocean  design  field.  Though  the  two  spectral  models  may  appear  very  different 
at  first  glance,  most  spectral  formulations  share  the  basic  form  of  Equation  2.6  [9] 

S(u>)  =  ^ H2se~B /"4  (2.6) 
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where 


A,  B  are  constants  to  be  determined 

Hs  is  the  significant  wave  height,  which  is  defined  to  be  the  average  wave  height  of  the 
one  third  highest  waves.  Sea  severity  is  most  commonly  given  in  terms  of  significant  wave 
height  among  researchers;  however,  significant  data  exists  which  gives  sea  severity  in  terms  of 
a  value  called  sea  state.  Appendix  B  provides  a  correlation  amongst  significant  wave  height, 
modal  wave  frequency,  wind  speed,  and  sea  state.  The  two  actual  derived  spectral  formulations 
chosen  for  the  model  are  the  Bretschneider  two  parameter  spectrum  and  the  Ochi  six  parameter 
spectrum,  given  in  Equations  2.7  and  2.8,  respectively.  [7] 

S(u)  =  h^.<AH (2.7) 


(^i±V 


mj/ 


H 2. 

S3 


r(A  j) 


U) 


4A.-+1 


=>-(— r~  )(w« 


iM4 


(2.8) 


The  other  two  models  are  the  square  wave  spectrum  and  the  single  sine  wave  spectrum, 
given  in  Equations  2.9  and  2.10,  respectively.  These  simple  geometric  form  models  were  used 
for  troubleshooting  and  debugging  of  the  model  and  to  verify  the  results  of  the  model. 


S(b)  = 


0.0  for  b  <  0.01 

5.0  for  0.01  <b  <1.01 
0.0  for  b  >  1.01 


(2.9) 


S  (6)  =  5 8{b-  bs) 


(2.10) 


where 
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Bretschneider  Energy  Density  Spectrum 


(De\elopment  =1.0  for  all  cases) 


Figure  2-3:  Bretschneider  Spectrum  for  Varying  Sea  Severity 


b=l  <»»> 

and  bs  is  the  spatial  frequency  of  the  sine  wave. 

2.2.1  Bretschneider  Two  Parameter  Spectrum 

The  two  parameters  of  the  Bretschneider  spectrum,  Equation  2.7,  are  significant  wave  height, 
Hs,  and  modal  frequency,  u>m.  Modal  frequency  is  the  frequency  at  which  S(u)  reaches  its 
maximum  value,  and  significant  wave  height  is  as  previously  defined.  For  the  special  case  of 
a  fully  developed  sea,  a  relationship  exists  between  um  and  Hs.  A  fully  developed  sea,  also 
referred  to  as  a  fully  arisen  sea,  is  one  in  which  the  wind  imparts  energy  to  the  waves  at  a 
rate  equal  to  that  at  which  viscous  damping  dissipates  their  energy  and  wave-wave  interactions 
transmit  energy  from  one  frequency  to  another.  Designating  this  modal  frequency  as  u>mFA ,  we 
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Bretschneider  Energy  Density  Spectrum 
(Hs  =  3.0  m  for  all  cases) 


Figure  2-4:  Bretschneider  Spectrum  for  Varying  Degrees  of  Development 


have 


WrnFA 


(2.12) 


which  establishes  the  energy  equilibrium  characteristic  of  a  fully  arisen  sea.  The  two  require¬ 
ments  for  a  fully  developed  sea  are 


•  sufficient  distance  over  which  the  wind  is  blowing  over  the  water,  or  fetch 

•  sufficient  length  of  time  the  wind  has  been  blowing,  or  duration. 

If  the  fetch  or  duration  are  insufficient,  the  modal  frequency  will  be  greater  than  that  given 
by  Equation  2.12;  if  the  winds  are  subsiding,  the  modal  frequency  will  be  less  than  that  given 
by  Equation  2.12.  This  is  best  understood  by  considering  the  mechanism  by  which  waves  are 
generated.  The  wind  causes  ripples,  or  capillary  waves,  on  the  sea  surface,  which  is  a  short 
wavelength,  high  frequency  phenomenon.  Thus  as  a  storm  is  building,  wave  energy  input  occurs 
at  high  frequency,  and  the  energy  gradually  is  transferred  through  wave-wave  interactions  to 
lower  frequencies.  Likewise,  as  a  storm  is  subsiding,  the  high  frequency  waves  attenuate  most 
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rapidly  due  to  viscous  effects,  and  the  result  is  a  skewing  of  the  spectral  energy  toward  lower 
frequencies.  By  introducing  a  term  we  will  call  the  development  factor,  Dev,  we  can  quantify 
the  degree  of  sea  arousal  as  given  in  Equation  2.13.  [9] 

U>m  =  Dev  ■  UJmFA  (2.13) 

Using  this  convention,  values  of  Dev  <1.0  correspond  to  a  decaying  sea,  and  values  of  Dev  >1.0 
correspond  to  a  building  sea.  Figures  2-3  and  2-4  show  representative  plots  of  S(u>)  as  significant 
wave  height  and  development,  respectively,  are  varied.  Note  that  in  Figure  2-4,  that  while 
changing  the  development  shifts  the  location  of  peak  spectral  energy  density,  it  keeps  the  area 
under  the  curve, and  thus  the  total  energy,  constant. 

2.2.2  Ochi  Six  Parameter  Spectrum 

The  Ochi  six  parameter  spectrum  is  a  generalization  of  the  Bretschneider  two  parameter  spec¬ 
trum  in  two  respects: 

•  two  spectral  peaks  are  present,  each  with  their  own  modal  frequency  and  significant  wave 
height,  and  their  contributions  add  together 

•  a  shape  factor,  A,  can  be  applied  to  each  spectral  peak  to  skew  the  curve  right  or  left, 
independent  of  the  significant  wave  height  or  modal  frequency.  [9] 

In  the  form  given  in  Equation  2.8,  each  parameter  is  an  independent  variable.  However, 
by  analyzing  over  800  wave  spectra  from  the  North  Atlantic,  Ochi  obtained  a  relationship  for 
the  six  parameters  in  terms  of  the  measured  significant  wave  height,  Hs.  [7]  The  result  of  this 
analysis  is  a  family  of  spectra  having  a  95%  confidence;  however,  the  results  given  in  Table 

2.2  represent  a  single  “most  probable”  spectrum  [7].  The  model  implements  this  single  most 
probable  form  for  obtaining  the  six  parameters  of  the  Ochi  spectrum.  A  sample  plot  of  the 
Ochi  most  probable  spectrum,  showing  its  characteristic  bimodal  shape,  is  shown  in  Figure  2-5. 
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(oas-^uu)  (<o)s 


Parameter 

Value 

Hs  i 

0.8AHS 

Hs2 

0MHs 

^ml 

0.70e~OM6Hs 

^m2 

l.l5e~om9H° 

Ai 

3.00 

1.54e-0062H> 

Table  2.2:  Ochi  Six  Parameter  Spectrum:  Most  Probable  Values 


Ochi  Energy  Density  Spectrum 
(Single  Most  Probable  Spectrum  for  Hg  =  3.0  m) 


Figure  2-5:  Ochi  Six  Parameter  Single  Most  Probable  Spectrum 
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Square  Wave  Energy  Density  Spectrum 


Spatial  Frequency  (m_1) 


Figure  2-6:  Square  Wave  Spectrum 

2.2.3  Square  Wave  Spectrum 

The  square  wave  spectrum  was  added  to  the  model  to  provide  for  verification  of  the  resultant 
multidirectional  wave  field  spectral  analysis.  By  providing  a  more  easily  recognizable  input 
spectrum,  the  goal  was  to  generate  a  more  easily  recognizable  output  spectrum.  A  plot  of  the 
square  wave  spectrum  is  shown  in  Figure  2-6. 

2.2.4  Single  Sine  Wave  Spectrum 

The  single  sine  wave  spectrum  was  added  to  the  model  to  provide  for  an  efficient  recognizable 
input  spectrum  which  could  be  used  for  troubleshooting  the  frequency  related  components  of 
the  multidirectional  wave  field.  By  providing  a  more  easily  recognizable  input  spectrum,  the 
goal  was  to  generate  a  more  easily  recognizable  output  spectrum.  A  plot  of  the  single  sine 
wave  spectrum  is  shown  in  Figure  2-7. 
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Single  Sine  Wa\e  Energy  Density  Spectrum 
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Spatial  Frequency  (m'1) 

Figure  2-7:  Single  Sine  Wave  Spectrum 

2.3  Unidirectional  Modeling  Methods 

As  mentioned  earlier,  the  wave  field  in  the  previous  model  consisted  of  a  unidirectional  random 
seaway  with  the  choice  of  only  two  tow  directions.  The  wave  field  in  the  current  model  makes 
use  of  the  unidirectional  method  to  construct  the  multidirectional  wave  field,  so  it  is  necessary 
to  completely  understand  the  unidirectional  method  developed  by  Ulrich  before  moving  on 
to  the  development  of  the  multidirectional  method.  The  discussion  and  derivation  of  the 
unidirectional  method  that  follows  comes  from  Ulrich.  [9] 

The  sea  domain  variables  of  interest  for  the  unidirectional  wave  field  are 

•  sea  elevation,  h 

•  fluid  horizontal  velocity,  u 

•  fluid  vertical  velocity,  w. 

Two  methods  were  employed  in  arriving  at  the  sea  domain  properties  of  the  unidirectional 
wave  field,  each  with  certain  advantages  and  disadvantages.  The  first  method  will  be  called 
the  Sine  Wave  Superposition  (SWS)  method,  and  the  second  method  will  be  called  the  Fast 
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Fourier  Transform  (FFT)  method.  The  source  code  for  both  of  these  computer  programs  can 
be  found  in  reference  [10]. 


2.3.1  Sine  Wave  Superposition  (SWS)  Method 

As  previously  mentioned  in  Equation  2.1,  the  sea  elevation  can  be  obtained  by  a  superposition 
of  sine  waves  having  a  random  phase  and  an  amplitude  prescribed  by  Equation  2.2.  In  like 
manner,  the  fluid  velocities  can  be  calculated  by  the  superposition  method.  The  formulas  for 
h,  u  and  w  in  an  earth-fixed  reference  frame  are  provided  in  Equations.  2.14  -  2.16. 


N 


h(xf,t) 

=  ^  (ii(io)  cos(ki(u)  x1  —  Uit  +  (ffj) 
i=i 

N 

(2.14) 

u{xf , t ) 

=  J2  ai(u)  u cos (kj(u>)  x'  —  Ujt  +  <ftj) 

i—1 

N 

(2.15) 

w(x\t) 

=  lj  sin(ki(u ))  xf  —  Uit  +  ^ ) 

1 

(2.16) 

By  invoking  a  coordinate  transformation  and  the  dispersion  relation  from  Table  2.1,  one  arrives 
at  the  forms  appropriate  to  the  coordinate  system  moving  with  the  antenna.  These  are  presented 
in  Equations  2.17  -  2.19. 


h(x,t) 


u(x,t ) 


w(x,t) 


N 


I  COS 


5>m 

i= 1 
N 

i— 1 
N 

^^a*(u;)  a;  sin 


'  j  v  ,  u>n 7.  ^ 

x  -  [uji - —  )t  +  (pi 


a;  cos 


i= i 


uj^U 

kJu)  x-(ui - 1 — )t  +  4>a 

9 

(jJ^U 

kAoo)  x  -  (u>i - - — )t  +  <& 

9 


(2.17) 

(2.18) 
(2.19) 


The  spectral  energy  density  function  S(u> )  used  in  this  formulation  was  a  single  sided  spectrum. 
The  advantages  of  the  SWS  method  are 


•  simplicity  in  coding  fluid  property  subroutines 

•  good  resolution  for  S(u>)  at  lower  frequencies,  where  most  of  the  wave  energy  resides 
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Time- Frequency 

Space-Spatial  Fre- 

Domain 

quency  Domain 

f  =  l/T 

b=l/L 

u)  —  2i rf 

k  —  27 xb 

k  =  (2? rf)2/g 

u>  =  \/2tx  gb 

Table  2.3:  Comparison  of  Fourier  Transform  Domain  Relationships 

•  frequency  meshing  (Aw)  is  unrelated  to  the  spatial  meshing  (Arc),  allowing  both  antenna 
element  length  and  wave  energy  frequency  resolution  to  be  arbitrarily  chosen. 


2.3.2  Fast  Fourier  Transform  (FFT)  Method 

This  method  takes  advantage  of  the  fact  that  Equations  2.14  -  2.16  are  a  Fourier  series,  and 
thus  the  functions  h,  n,and  w  are  the  inverse  Fourier  transform  of  some  expression.  Since  the 
calculations  to  determine  the  sea  domain  properties  occur  at  a  fixed  point  in  time,  the  transform 
variables  are  not  u  and  t,  rather  they  become  b  and  x.  Table  2.3  summarizes  the  analogous 
relationships  between  the  time-frequency  domains  and  the  space-spatial  frequency  domains. 
Because  the  FFT  algorithm  used  for  the  model  assumed  the  presence  of  both  positive  and 
negative  frequency  components,  it  was  necessary  to  use  a  two  sided  spectrum  for  this  method. 
In  order  to  retain  the  same  amount  of  energy  (i.e.  area  under  the  spectral  energy  density  curve) 
as  the  transition  from  one-sided  to  two-sided  spectra  is  made,  it  is  necessary  that  the  formula 
for  the  component  amplitudes  be  adjusted  from  that  given  in  Equation  2.2.  In  addition,  it  is 
convenient  to  incorporate  the  random  phase,  <f>i}  into  a  complex  component  amplitude,  c*.  With 
these  changes  in  mind,  the  formulas  providing  the  sea  domain  properties  are  given  in  Equations 
2.20  -  2.24. 


h(x )  =  Re 
u(x)  =  Re 
w(x)  =  Re 


N  _  'I 

yi  Cj(bj)  exp  ji(27r6iX  +  {2ixbiU  -  y/2irgbi)  t)l  > 

i=l  J 

{N 

^  <k(bi)y/2ngbi  exp  \i{2'Kbix  +  (27T 6*17-  y/2i xgbi)t) 

i=i  ) 

N  _  _ 

y:  -i  Ci(bi)  y/lngbi  exp  ^(2^6^  +  (2nbiU  -  a/2tt gbi)  t) 


} 


(2.20) 

(2.21) 

(2.22) 
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Re{ci(6j)}  = 


is(£0A6  008(0,) 


(2.23) 


lm{ci(bi)}  =  S{bi)Ab  sin(0,) 
Then  the  Fourier  transform  pairs  become 

H(b)  &  h{x ) 
\Z2ngbi  H(b)  &  u(x) 
—iyj2'i\gbl  H(b)  w(x) 


where 

N 

H(b)  =  ^c,(h)exp  j it  (2irbiU  -  y/2ngb^ 

i= 1 

and 


(2.24) 


(2.25) 

(2.26) 
(2.27) 


(2.28) 


H(b)  &  h(x) 

is  defined  to  mean 

N 

m  =  £m  x)  exp  {—i2-Kbx) 
i= 1 

1  N 

h(x)  =  —  H (b)  exp  ( i2irbx ) 

2?r  i=i 

By  setting  the  complex  component  amplitudes,  c,,  to  obey  the  relation 


Ciipi)  —  Cj(  ^i) 


(2.29) 

(2.30) 

(2.31) 

(2.32) 


that  is,  the  positive  frequency  components  are  the  complex  conjugates  of  the  negative  frequency 
components,  then  the  imaginary  parts  of  Equations  2.20-2.22  are  zero,  and  the  Re{  }  can  be 
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dropped  from  the  equations.  Finally,  the  wave  field  is  constructed  by  performing  an  inverse 
FFT  on  the  one-dimensional  complex  arrays,  represented  by  the  left  side  of  Equations  2.25  - 
2.27,  which  contain  the  Fourier  coefficients  describing  h,  u,  and  w. 

When  converting  between  the  spectral  energy  density  function  as  a  function  of  u  to  another 
variable,  such  as  b,  the  basic  relations  of  Equations  2.33  and  2.34  are  useful. 

S{u)dw  =  S(b)  db  (2.33) 

S(h)  =  S(  u)~  (2.34) 

An  important  feature  of  the  FFT  method  is  the  linkage  of  the  spatial  frequency  domain 
meshing  (A6),  and  the  sea  spatial  meshing  (Ax).  The  key  to  the  FFT  method  is  to  make  the 
antenna  spatial  meshing,  As,  match  the  sea  domain  spatial  meshing,  Ax.  This  is  accomplished 
by  using  the  relations  given  in  Equations  2.35  -  2.38,  and  ensuring  that  the  sea  Ax  equals  the 
desired  antenna  As. 


A  b  = 

i  , 

—  °min 

^sea 

(2.35) 

frmax  — 

NSea^b  Ngea 

l 

(2.36) 

2  2  Lsea 

-  2  Ax 

^min  ~ 

^/Zngbmm 

(2.37) 

^max  = 

^/Zngbmzx 

(2.38) 

The  primary  advantage  of  the  FFT  method  is  computational  speed.  The  FFT  method 
executes  approximately  10  times  faster  than  the  SWS  version.  A  further  increase  in  execution 
speed  is  accomplished  by  not  computing  the  FFT  every  time  step,  At.  Rather,  the  sea  data 
(h,  u,  w)  is  computed  at  nAt  intervals,  and  linearly  interpolated  in  between.  The  integer  n  is 
determined  by  the  computer  model  using  the  highest  sea  energy  frequency  component,  wmax. 
This  ensures  that  the  linearly  interpolated  sea  evolves  in  a  smooth  and  continuous  fashion. 

Its  disadvantage  is  that  it  has  poorer  low  frequency  resolution  of  the  sea  wave  energy  than  its 
SWS  counterpart.  This  is  demonstrated  by  a  typical  run  case  in  Table  2.4  where  the  lowest  10 
frequencies  defined  are  listed  for  both  methods.  As  one  can  readily  calculate  from  integrating 
Equation  2.7  with  respect  to  ui,  for  a  sea  severity  of  2.0  m  significant  wave  height,  67.7%  of  the 
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Component 

Frequency,  w, 
SWS  Method 

FFT  Method 

i 

0.347 

0.347 

2 

0.381 

0.491 

3 

0.414 

0.601 

4 

0.448 

0.694 

5 

0.482 

0.776 

6 

0.515 

0.850 

7 

0.549 

0.918 

8 

0.583 

0.981 

9 

0.617 

1.041 

10 

0.650 

1.097 

Notes: 

1.  FFT  case: 
components 

to  €  [0.347,22.2]  with  4096  frequency 

2.  SWS  case: 
components 

w  G  [0.347,24.0]  with  700  frequency 

Table  2.4:  Comparison  of  Frequency  Resolution  for  the  SWS  and  FFT  Methods 


sea  wave  energy  lies  below  a  frequency  of  1.097  rad/sec.  Thus  the  FFT  method  suffers  in  that 
it  places  most  of  the  frequency  components  at  high  frequencies,  where  little  of  the  wave  energy 
resides.  This  undesirable  result  is  a  consequence  of  choosing  a  b  space  meshing  with  constant 
A 6,  which  is  necessary  to  compute  FFT’s  in  the  b  and  x  domains.  A  constant  A b  meshing 
produces  a  non-constant  Aw  meshing  with  Aw*  =  \/‘2ir()bn]\u  (Vi  —  V%~  l)  • 


2.4  Multidirectional  Modeling  Methods 

The  current  model  improves  upon  the  previous  model  by  upgrading  the  two-dimensional  unidi¬ 
rectional  wave  field  discussed  in  Section  2.3  to  a  three-dimensional  multidirectional  wave  field 
with  the  choice  of  any  arbitrary  tow  angle  in  the  wave  field. 

The  sea  domain  variables  of  interest  for  the  multidirectional  wave  field  are 

•  sea  elevation,  h 

•  horizontal  fluid  velocity  (x-direction) ,  u 

•  horizontal  fluid  velocity  (y-direction) ,  v 
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•  vertical  fluid  velocity,  w. 


Two  methods  were  considered  in  arriving  at  the  sea  domain  properties  of  the  multidirectional 
wave  field,  each  with  certain  advantages  and  disadvantages.  The  first  method  will  be  called 
the  Complete  Wave  Field  (CWF)  method,  and  the  second  will  be  called  the  Partial  Wave  Field 
(PWF)  method. 

2.4.1  Complete  Wave  Field  (CWF)  Method 

The  Complete  Wave  Field  (CWF)  method  involves  the  construction  of  a  three-dimensional  wave 
field  using  a  spectral  energy  model  that  is  completely  described  in  both  the  x  and  y  directions 
in  the  spatial  frequency  domain.  This  is  accomplished  in  a  manner  similar  to  that  discussed 
in  Section  2.3.2,  except  that  the  equations  now  will  be  a  function  of  both  x  and  y  and  will 
separate  the  horizontal  velocity  into  x  and  y  components. 

As  previously  mentioned  in  Equation  2.1,  the  sea  elevation  can  be  obtained  by  a  superpo¬ 
sition  of  sine  waves  having  a  random  phase  and  an  amplitude  prescribed  by  Equation  2.2.  In 
like  manner,  the  fluid  velocities  can  be  calculated  by  the  superposition  method.  The  formulas 
for  h,  u,  v,  and  w  in  an  earth-fixed  reference  frame  are  provided  in  Equations  2.39  -  2.42. 

N  N 

h{x',  y\  t)  =  aJk(uXj , uVk)  cos (kXj (uXj )  x'  +  kyk  (uyk )  y'  -  ojjkt  +  <pjk)  (2.39) 

3= 1  k= 1 

N  N 

u(x',y',t)  =  XI  aik(uXj ,  0Jyk)  oJjk  cos  {6jk)  (2.40) 

j= l  fc=i 

x  cos {kXj (ojXj )  x'  +  kyk  (u jyk)  y'  -  ujkt  +  <f>jk) 

N  N 

v(x',y',t)  =  ajk(oJXj ,  0Jyk )  ujk  sin  (9jk)  (2.41) 

j= 1  k= 1 

x  cos {kXj (uXj)x'  -|-  kVk  (ujyk)  y'  -  ojjkt  +  (f>jk) 
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N  N 

w(x',y',t)  =  Y^Yiajk^Xj,uyk)  Ujksm{kXj{uXj)x'  +  kyk(u)yk)y' -Ujkt  +  <!>jk)  (2.42) 

j= l  fc=i 

where 


6jk  =  tan  1  f 

\fcxj  J 


(2.43) 


(2.44) 


By  invoking  a  coordinate  transformation  and  the  dispersion  relation  from  Table  2.1,  one 
arrives  at  the  forms  appropriate  to  the  coordinate  system  moving  with  the  antenna.  These  are 
presented  in  Equations  2.45  -  2.48. 


h(x,y,t) 


N  N 


j- 1  fc=l 


(2.45) 


X  cos 


u{x,y,t) 


N  N 

}  ^  ^  ajkiyixj  5  Myk)  ^ jk  cos  (djk) 

j= 1 


(2.46) 


X  cos 


uj^U  cos  (djk) 

kxji^Xj)  x  +  kyk  (uiyk)  y  —  {bJjk  ~  )t  +  </) jk 


v{x,y,t) 


N  N 

ajk{^xj  >  Wyk )  ujjk  sin  (djk) 

j=l  k=  1 


(2.47) 


X  cos 


u%U  cos  ( Ojk ) 

kXj (ujXj )x  +  kyk  (uyk )  y  -  (ujk - - - )t  +  (j>jh 
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Bretschneider  Energy  Density  Spectrum 


Figure  2-8:  Complete  Wave  Field  3-D  Energy  Spectrum  (quadrants  1  and  4  shown) 


w(x,y,t) 


N  N 


j= i  k=i 


x  sin 


+  kyk  (u  Vk  )  y  -  (“jk 


(2.48) 


uj2jkU  cos  {6jk) 
9 


)t  +  fijk 


We  again  take  advantage  of  the  fact  that  Equations  2.45  -  2.48  are  a  Fourier  series,  and 
thus  the  functions  h,  u,  u,and  w  are  the  inverse  Fourier  transform  of  some  expression.  Since 
the  calculations  to  determine  the  sea  domain  properties  occur  at  a  fixed  point  in  time,  the 
transform  variables  are  not  u  and  t ,  rather  they  become  b  and  x.  Because  the  two-dimensional 
FFT  algorithm  used  for  the  model  assumed  the  presence  of  frequency  components  in  all  four 
quadrants  surrounding  the  origin  in  the  spatial  frequency  domain,  it  was  necessary  to  construct 
an  energy  density  spectrum  in  three  dimensions.  An  example  of  a  three-dimensional  energy 
density  spectrum,  in  spatial  frequency  coordinates,  which  is  used  by  the  sea  domain  properties 
equations  is  shown  in  Figure  2-8.  In  order  to  retain  the  same  amount  of  energy  (i.e.  area  under 
the  spectral  energy  density  curve),  it  is  necessary  that  the  formula  for  the  component  amplitudes 
be  adjusted  from  that  given  in  Equation  2.2.  In  addition,  it  is  convenient  to  incorporate  the 
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random  phase,  into  a  complex  component  amplitude,  Cjk ■  With  these  changes  in  mind,  the 
formulas  providing  the  sea  domain  properties  are  given  in  Equations  2.49  -  2.54. 


h(x,y) 


N  N 


y '  y !  Re  { cjk{bXj  >  byk ) 

j=i  /:= i 


(2.49) 


x  exp 


i(2nbXjx  +  271 \ky  +  {2irbjkU cos  (0jk)  -  y/2irgbjk)  t )  j 


N  N  _ 

u(x,y)  =  X  X Re  [cjkibx^by^^ngbjk  cos  (0jk)  (2.50) 

j= l  fc= l 

x  exp  ^i(2irbXjx  4-  2-xbyky  +  (2nbjkU cos  (0jk)  -  ^J2'ngb]k)  i)]  } 


N  N 

v(x,y )  =  X  X  Re  {^(bxj ,  byk)^2ngb~k  sin  (0jk)  (2.51) 

j= l  fc= l 

x  exp  [i(2TrbXjx  +  2nbVky  +  {2nbjkU  cos  (0jk)  -  y/2ngbjk)  t)J  } 


w(x,y) 


N  N 

X  X  Re  {  crt(bxj  >  byk )  \/27 Tgb~k 

j=  1  k= 1 


(2.52) 


x  exp  i(2nbXjx  +  2rcbVky  +  ( 2-nbjkU  cos  ( 6jk )  -  ^2-ngbjk)  t )  | 


R e{cjk(bXj,byk)}  =  ^5(6,,^  Abcos(<J>jk) 
lm{cjk(bXj,byk)}  =  ]j^S(bXj,~b^)  Absm((f>jk) 


(2.53) 


(2.54) 


where 


h*  =  +  % 


(2.55) 
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Then  the  Fourier  transform  pairs  become 


H  ( bx ,  by) 

<=> 

h{x,  y) 

(2.56) 

\/2i rgbjk  cos  (9jk)  H(bx ,  by) 

u{x,y) 

(2.57) 

y/2i xgbjk  sin  (0jk)  H(bx,by) 

v(x,y) 

(2.58) 

—i-\](2'Kgbjk  H  (bx,  by) 

w(x,y) 

(2.59) 

where 


N  N 


H(bx,  by)  =  cjk(bx^byk)  exp  it  (2nbjkU  cos  (9jk)  -  y/2irgbjkj  (2.60) 


3=1  fc=i 


and 


H(bx,by)  h(x,y)  (2.61) 

is  defined  to  mean 

N  N 

H(bx ,  by)  =  X)  h(xj>yk)  exP  {-i2nbXjx  ~  i2irbyky)  (2.62) 

j=i  k= 1 
^  N  N 

h^y)  =  ^EE^*  byic )  exp  (i2nbXj x  +  i2nbyky)  (2.63) 

X  i=i  fc=i 

By  setting  the  complex  component  amplitudes,  Cjk,  to  obey  the  relations 


cjk{  bx. ,  6yfc)  —  Ci {pxj )  byk )  (2-64) 

Cjk{~bXj  ■  byk)  =  Ci{bxj>~byk )  (2.65) 

that  is,  the  frequency  components  of  quadrants  2  and  3  are  the  complex  conjugates  of  the 
frequency  components  of  quadrants  4  and  1,  respectively,  then  the  imaginary  parts  of  Equations 
2.49  -  2.52  are  zero,  and  the  Re{  }  can  be  dropped  from  the  equations. 

Finally,  the  wave  field  is  constructed  by  performing  an  inverse  FFT  on  the  two-dimensional 
complex  arrays,  represented  by  the  left  side  of  Equations  2.56  -  2.59,  which  contain  the  Fourier 
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Three  Dimensional  Wave  Field  -  CWF  Method 
NOTE:  Only  Portion  of  Wave  Field  Shown 


Y(meters)  0  °  X  (meters) 


Figure  2-9:  Wave  Field  Generated  by  CWF  Method 

coefficients  describing  and  w.  A  example  of  the  wave  field  generated  by  the  CWF 

method  is  shown  in  Figure  2-9.  This  computational  step  is  not  trivial  and  is  the  reason  that 
this  method  was  abandoned  early  in  the  project.  The  memory  requirements  and  computation 
times  necessary  to  complete  the  inverse  FFT  of  these  large  arrays  was  not  practical  for  desktop 
computer  modeling. 

Wave  field  directionality  has  not  been  discussed  up  to  this  point.  That  topic  has  been 
reserved  for  Section  2.4.2,  the  Partial  Wave  Field  method,  which  is  the  method  that  is  imple¬ 
mented  in  the  model. 

2,4.2  Partial  Wave  Field  (PWF)  Method 

The  Partial  Wave  Field  (PWF)  method  involves  the  construction  of  a  three-dimensional  wave 
field  using  a  spectral  energy  model  that  is  only  partially  described  in  both  the  x  and  y  directions 
in  the  spatial  frequency  domain.  In  comparison  to  the  CWF  method,  this  partial  description 
significantly  reduces  the  amount  of  computation  needed  to  build  the  modeled  wave  field.  It 
reduces  computation  time  by  forming  the  multidirectional  wave  field  from  the  proper  combi¬ 
nation  of  several  unidirectional  waves  as  shown  in  Figure  2-10.  In  addition,  the  PWF  method 
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Figure  2-10:  PWF  Unidirectional  Wave  Summation 

takes  advantage  of  the  fact  that  the  only  wave  field  information  necessary  to  compute  exposure 
statistics  on  the  antenna  is  the  information  from  the  wave  field  that  is  directly  in  contact  with 
the  antenna. 

The  first  step  of  the  PWF  method  is  to  select,  from  the  available  choices,  the  wave  energy 
density  spectrum  to  be  modeled.  The  energy  density  spectrum  dictates  the  behavior  and 
magnitude  of  the  modeled  seas  and  its  selection  can  significantly  affect  the  exposure  statistics. 
Guidance  on  proper  selection  of  the  energy  density  spectrum  can  be  found  in  Section  A.2.1  and 
in  numerous  texts  and  papers.  [4]  [5]  [7]  [8] 

The  next  step  of  the  PWF  method  is  to  construct  the  three  dimensional  energy  density 
spectrum  with  NN  discrete  unidirectional  spectra  that  are  each  propagating  in  a  different 
direction.  The  model  allows  for  the  user  to  select  any  odd  number  of  discrete  unidirectional 
waves.  The  requirement  for  an  odd  number  is  driven  by  the  assumption  that  one  direction  lies 
on  the  principal  axis  and  the  remaining  directions  are  split  evenly  on  either  side  of  the  principal 
axis.  A  larger  number  of  directions  will  result  in  a  more  realistic  three  dimensional  sea  spectrum 
but  at  the  expense  of  computation  time.  Referring  again  to  Figure  2-10,  one  can  see  that  as 
the  number  of  unidirectional  waves  increases,  the  modeled  seas  more  closely  approach  a  realistic 
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Bretschneider  Energy  Density  Spectrum 


Figure  2-11:  Partial  Wave  Field  3-D  Energy  Spectrum  -  No  Directional  Scaling  (quadrants  1 
and  4  shown) 

three  dimensional  spectrum.  Likewise,  as  the  number  of  unidirectional  waves  decreases,  the 
modeled  seas  move  further  away  from  a  realistic  three  dimensional  spectrum.  The  model 
constructs  the  three  dimensional  energy  density  spectrum  by  assigning  a  direction  to  each  of 
the  discrete  unidirectional  spectra,  as  shown  in  Figure  2-11.  The  example  shown  in  Figure 
2-11  is  for  five  propagation  directions.  For  a  given  number  of  propagation  directions,  five  in 
this  example,  the  model  assumes  equal  angles  between  each  direction  and  then  solves  for  the 
correct  angle  to  assign  to  each  direction.  It  should  be  noted  that,  for  purposes  of  illustration, 
the  magnitude  of  each  of  the  unidirectional  spectra  is  equal  and  that  Figure  2-11  is  similar 
to  Figure  2-8  except  that  Figure  2-11  is  only  partially  modeling  the  three  dimensional  energy 
density  spectrum. 

It  has  already  been  pointed  out  the  magnitude  of  each  of  the  unidirectional  spectra  in  Figure 
2-11  is  equal.  This  implies  that  the  wave  field  has  no  directional  spreading  characteristic;  that 
is,  energy  is  being  propagated  equally  in  all  directions.  In  reality,  wind-generated  waves  tend 
to  propagate  in  the  same  general  direction  as  the  wind.  Typically,  the  energy  is  spread  over 
several  directions,  but  the  majority  is  in  the  direction  of  the  wind.  This  spreading  of  energy 
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can  be  characterized  by  what  is  known  as  a  directional  spreading  function.  The  directional 
spreading  function  assigns  a  fraction  of  the  total  wave  energy  to  each  particular  direction,  with 
that  fraction  being  determined  by  the  shape  of  the  spreading  model  selected.  To  generate 
a  directional  spectrum,  one  simply  multiplies  the  energy  density  spectrum  by  the  directional 
spreading  function 


S(v,0)  =  S(u>)D(Otu)  (2.66) 

where  D  ( 9 ,  u)  represents  the  directional  spreading  as  a  function  of  both  direction,  0 ,  and  fre¬ 
quency,  u>.  Not  unlike  wave  spectra,  numerous  directional  spreading  functions  have  been  pro¬ 
posed  by  researchers  attempting  to  characterize  a  particular  sea  wave  environment.  [4]  [5]  [6]  [7] [8] 
Although  the  spreading  functions  implemented  in  the  model  are  simple  models  that  are  a  func¬ 
tion  of  direction  only,  spreading  models  exist  that  represent  the  spreading  as  a  function  of  both 
direction  and  frequency.  In  these  frequency  dependent  spreading  functions,  the  low  frequency 
components  are  assigned  a  very  narrow  spreading  while  the  high  frequency  components  are  as¬ 
signed  a  wider  spreading.  [5]  The  spreading  functions  used  in  the  model  are  the  cosine  square 
formula  and  the  cosine  fourth  formula,  given  in  Equations  2.67  and  2.68,  respectively. 

D(9)  =  ^cos2(0),  -1<0<1  (2.67) 

D{6)  =  £«*(*),  -1<9<!  (2.68) 

A  plot  of  the  cosine  square  formula,  shown  in  Figure  2-12,  shows  the  behavior  of  the  di¬ 
rectional  spreading  function  as  the  propagation  direction  changes  from  0  degrees,  the  principal 
axis  of  propagation,  to  ±90  degrees.  As  shown  in  Figure  2-12,  the  largest  fraction  of  energy 
is  along  the  principal  axis  while  at  ±90  degrees  the  fraction  tapers  to  zero.  To  account  for 
the  total  energy  in  the  spreading  spectrum,  the  actual  fraction  assigned  to  each  unidirectional 
spectrum  is  the  integral  of  the  directional  spreading  function  centered  on  each  direction  with 
integration  limits  being  equal  distances  on  either  side  of  the  propagation  direction.  Again, 
using  the  five  propagation  direction  example,  the  propagation  directions  and  their  respective 
integration  limits  are  shown  in  Figure  2-13. 

Applying  the  cosine  square  spreading  function  to  the  three  dimensional  energy  density  spec- 
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(*/  ) cos20  Directional  Spreading  Function 


Figure  2-12:  Cosine  Square  Directional  Spreading  Function 


(2/  )cos20  Directional  Spreading  Function 


Figure  2-13:  Cosine  Square  Directional  Spreading  Function  for  the  Five  Directions  Example 
Problem 


38 


Bretschneider  Energy  Density  Spectrum 


Figure  2-14:  Partial  Wave  Field  3-D  Energy  Spectrum  -  With  Directional  Scaling  (quadrants  1 
and  4  shown) 

trum  in  Figure  2-11,  you  obtain  a  three  dimensional  energy  density  spectrum  with  directional 
characteristics,  shown  in  Figure  2-14.  Now,  comparing  Figure  2-14  with  Figure  2-11,  you  can 
see  that  the  magnitudes  of  the  unidirectional  spectra  in  Figure  2-14  are  not  equal,  but  are  scaled 
according  to  the  directional  spreading  function.  In  addition,  as  prescribed  by  the  spreading 
function,  the  principal  axis  of  wave  propagation,  0  degrees,  contains  the  largest  fraction  of 
energy. 

The  next  step  in  the  PWF  method  is  to  apply  the  FFT  method,  described  in  Section  2.3.2, 
to  each  of  the  scaled  unidirectional  wave  spectra  to  solve  for  the  spatial  sea  domain  properties 
of  each  unidirectional  wave.  Now  the  individual  waves  must  be  properly  summed  to  determine 
the  three  dimensional  wave  properties  along  the  antenna.  The  sea  domain  properties  of  the 
PWF  method  are  somewhat  different  than  those  found  by  the  CWF  method.  In  the  CWF 
method,  the  goal  was  to  solve  for  h,  u,  v,  and  w  in  the  entire  wave  field.  With  the  PWF 
method,  the  only  information  needed  is  the  information  along  the  antenna.  To  simplify  the 
numerical  modeling,  the  velocities  of  interest  were  found  relative  to  the  antenna.  So  rather 
than  solving  for  the  x  and  y  velocities  in  the  entire  field,  the  model  solves  for  the  velocities 
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that  act  longitudinally  and  transversely  on  the  antenna,  as  well  as  the  vertical  velocity.  The 
summation  is  accomplished  by  Equations  2.69  -  2.72. 

NN 

h  (x)  =  ^2  h  (£)  (2.69) 

1=1 


NN 

^ longitudinal  (^)  =  £  Ui  ( x )  cos  (ip)  (2.70) 

2=1 


NN 

M transverse  (®)  =  £  Ui  (x)  sin  (ip)  (2.71) 

2  =  1 


NN 

w  (x)  =  5^  Wi  (®)  (2.72) 

i= 1 

where 


x  —  -  0.5 1  As  -  -  [|(JV-l)  +  l]  +0.5  j  As  cos  (VO  (2.73) 

and 

NN  is  the  number  of  propagation  directions  being  summed 

hi  ( x )  is  the  wave  elevation  of  the  individual  unidirectional  wave,  at  the  location  x,  in  the 
unidirectional  wave  field  that  is  acting  on  the  antenna  location  x,  as  shown  in  Figure  2-15 
Ui  (x)  is  the  horizontal  velocity  of  the  individual  unidirectional  wave,  at  the  location  x,  in 
the  unidirectional  wave  field  that  is  acting  on  the  antenna  location  x,  as  shown  in  Figure  2-15 
ip  is  the  angle  between  the  individual  unidirectional  wave  and  the  antenna,  as  shown  in 
Figure  2-15 

Wi  (x)  is  the  vertical  velocity  of  the  individual  unidirectional  wave,  at  the  location  x,  in  the 
unidirectional  wave  field  that  is  acting  on  the  antenna  location  x,  as  shown  in  Figure  2-15 
Nsea  is  the  number  of  spatial  divisions  in  the  unidirectional  wave  field 
As  is  the  spatial  increment  on  the  antenna  as  well  as  the  spatial  frequency  increment  of  the 
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Figure  2-15:  Geometry  of  the  Wave  Summation 

unidirectional  wave  spectrum 

N  is  the  number  of  spatial  divisions  on  the  antenna 
L  is  the  modeled  length  of  the  antenna 

In  Equations  2.69  -  2.73,  it  is  assumed  that  the  width  of  the  unidirectional  wave  is  infinite 
and  that  the  surface  properties  along  that  infinite  width  do  not  vary.  This  assumption  allows 
you  to  perform  the  summation  shown  in  Figure  2-15. 

The  origin  of  equation  2.73  is  not  obvious,  but  it  represents  the  trigonometric  and  geometric 
relation  between  the  individual  wave  and  the  antenna.  Referring  to  Figure  2-15,  the  numbers 
represent  the  array  storage  locations  of  the  particular  variables  with  the  wave  properties  num¬ 
bered  i  =  1  to  Nsea  and  the  antenna  properties  numbered  j  —  1  to  N.  The  spacing  between 
points,  As,  is  the  same  for  both  the  antenna  and  the  wave  with  the  length  of  the  wave  field 
being  Lsea  and  the  length  of  the  antenna  being  L.  In  addition,  it  is  assumed  that  the  antenna 
and  the  wave  cross  at  their  midpoints.  The  position,  x ,  can  be  found  from 
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where 


^crossing  point 


^origin 


=  (number  of  whole  and  partial  mesh  spaces  on  wave  field  between 
crossing  point  and  origin)  x  (As) 

=  jp  -  0.5^  As  (2.75) 

=  0.0  (2.76) 


^crossing  point 


X 


(number  of  whole  and  partial  mesh  spaces  on  antenna  between 
crossing  point  and  origin)  x  (As) 

(f-  0.5^  As  (2.77) 

(number  of  whole  and  partial  mesh  spaces  on  antenna  between 
x  and  origin)  x  (As)  (2.78) 

0-1)  As  (2.79) 


combining  equations  2.74  -  2.79  you  obtain 

N sea 


X  = 


-  0.5^  As  -  -  j  +  0.5^  As  cos  (ip) 


(2.80) 


where  j  represents  the  array  storage  location  nearest  to  x,  given  by 


i  =  |(j\r-i)  +  i 


(2.81) 


when  Equation  2.81  is  combined  with  Equation  2.80,  you  obtain 


LZ(iV-l)  +  l 


]+0.5} 


As  cos  (ip) 


(2.82) 


which  is  identical  to  Equation  2.73. 

The  PWF  modeling  framework  consists  of  the  following  steps 

1.  select  the  wave  energy  density  spectrum  to  be  modeled 

2.  construct  the  three  dimensional  energy  density  spectrum  with  NN  discrete  unidirectional 
spectra  that  are  each  propagating  in  a  different  direction 
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3.  select  the  directional  spreading  function 

4.  apply  the  directional  spreading  function  to  the  discrete  components  of  the  three  dimen¬ 
sional  energy  density  spectrum 

5.  perform  an  inverse  FFT  on  the  one-dimensional  complex  arrays,  represented  by  the  left 
side  of  Equations  2.25  -  2.27,  which  contain  the  Fourier  coefficients  describing  h,u,  and 
w  of  each  discrete  unidirectional  spectrum 

6.  properly  sum  the  sea  domain  properties  of  each  discrete  unidirectional  wave  along  the 
length  of  the  antenna  to  generate  a  three  dimensional  wave  field  on  the  antenna 

7.  time  advance  the  Fourier  coefficients  describing  h,  u,  and  w  using  Equation  2.28 

8.  repeat  steps  5-7  for  duration  of  simulation. 

2.5  Verification  of  Model  Seas 

The  random  multidirectional  seaway  was  constructed  by  the  PWF  method  as  described  in 
Section  2.4.2.  By  running  the  modeling  code  and  sampling  the  wave  field  for  several  parameters, 
one  can  build  a  wave  record  that  can  be  analyzed  to  determine  if  the  resultant  wave  field 
correctly  models  the  desired  wave  field.  The  modeled  seas  verification  focused  on  the  following 
four  areas 

•  significant  wave  height,  Hs 

•  encounter  frequency,  toe 

•  spectral  analysis  of  a  spatial  wave  record 

•  spectral  analysis  of  a  temporal  wave  record. 

2.5.1  Significant  Wave  Height 

By  sampling  the  wave  record  for  sea  surface  elevation,  one  can  determine  the  significant  wave 
height  of  the  sample,  as  defined  in  Section  2.2.  [7]  This  can  then  be  compared  to  the  significant 
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Sample 

Mean(/ij) 

(m) 

Spatial 

a 

(m) 

Hs 

(m) 

Temporal 

Mean(/ij)  a 

(m)  (m) 

Hs 

(m) 

i 

-0.0546 

0.7177 

2.8779 

0.0102 

0.7784 

3.1215 

2 

-0.0120 

0.7455 

2.9893 

0.0032 

0.6950 

2.7869 

3 

0.0024 

0.8389 

3.3638 

-0.0009 

0.7893 

3.1652 

4 

0.0160 

0.8890 

3.5647 

0.0085 

0.7582 

3.0405 

5 

0.0287 

0.6336 

2.5406 

0.0231 

0.7017 

2.8137 

6 

-0.0119 

0.8843 

3.5459 

-0.0045 

0.7159 

2.8707 

7 

0.0066 

0.7410 

2.9713 

0.0057 

0.7528 

3.0186 

8 

-0.0844 

0.7636 

3.0621 

-0.0171 

0.7140 

2.8632 

9 

-0.0362 

0.5490 

2.2016 

-0.0185 

0.7412 

2.9723 

10 

-0.0391 

0.8569 

3.4360 

0.0092 

0.7373 

2.9567 

11 

0.0632 

0.7047 

2.8258 

-0.0182 

0.7215 

2.8934 

12 

0.0050 

0.6758 

2.7100 

-0.0061 

0.7636 

3.0619 

Averages: 

-0.0097 

0.7500 

3.0075 

-0.0005 

0.7391 

2.9637 

Table  2.5:  Simulated  Significant  Wave  Heights  for  a  Specified  Sea  Severity  of  3.0  m 


wave  height  which  was  specified  in  the  input  file.  The  measured  significant  wave  height  can  be 
found  by  Equations  2.83  and  2.84, 


1  71 

cr  =  ,  >  x\ 

n  —  1 

7—1 

(2.83) 

Hs  =  4.01 

(2.84) 

where 

Xi  is  the  sampled  wave  amplitude,  measured  from  the  mean  free-surface 
a  is  the  standard  deviation  of  the  sampled  wave  amplitudes 
H$  is  the  measured  significant  wave  height  of  the  wave  record. 

(The  “  ^  ”  symbol  denotes  that  the  quantity  is  computed  from  a  sampled  wave  record,  and 
therefore  may  deviate  somewhat  from  the  theoretical  value.) 

A  sea  severity  of  3.0  m  significant  wave  height  was  simulated  by  the  PWF  method.  Twelve 
wave  amplitude  records  were  analyzed  by  sampling  the  seas  spatially  at  0.0625  rrt  intervals 
and  temporally  at  0.021  sec  intervals.  The  results  are  summarized  in  Table  2.5.  An  excellent 
agreement  to  the  specified  significant  wave  height  is  achieved  with  an  average  deviation  in  Hs 
of  0.25%  and  1.20%  for  the  spatial  and  temporal  analysis,  respectively. 
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Encounter  Frequency  for  Varying  Tow  Speeds 


Figure  2-16:  Encounter  Frequency  for  Varying  Tow  Speeds 


2.5.2  Encounter  Frequency 

Determination  of  the  encounter  frequency  presents  a  formidable  task  for  a  random  multidirec¬ 
tional  seaway.  To  simplify  the  task,  the  analysis  of  encounter  frequency  was  done  using  the 
single  sine  wave  spectrum  as  the  input  spectrum.  The  method  used  to  verify  the  encounter 
frequency  invoked  the  superposition  principle:  the  sea  waves  generated  by  the  PWF  method  are 
ultimately  a  superposition  of  individual  sine  waves,  thus  the  verification  of  encounter  frequency 
was  restricted  to  seas  represented  by  a  single  sine  wave.  The  goal  of  this  analysis  is  to  verify 
that  the  antenna  sees  the  correct  “encounter  frequency,”  i.e.  an  effective  change  of  the  wave 
frequency  as  the  antenna  is  towed  in  a  seaway.  The  encounter  frequency,  u>e,  is  given  by 


We  =  U  0 


ojqU  cos  (6) 
9 


(2.85) 


where 

wo  is  the  wave  frequency  in  a  stationary  coordinate  system 
U  is  the  antenna  tow  speed 

9  is  the  tow  angle  with  respect  to  the  principal  axis  of  wave  propagation 
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Spatial  Power  Spectral  Density  Analysis 


1  Propagation  Direction 


Figure  2-17:  Spectral  Analysis  of  Spatial  Wave  Data  -  One  Propagation  Direction 
g  is  the  acceleration  due  to  gravity. 

Thus  for  towing  into  head  seas  (cos  (180°)  is  negative),  the  encounter  frequency  goes  up, 
and  for  towing  with  following  seas  (cos(0°)  is  positive),  the  encounter  frequency  goes  down. 
The  encounter  frequency  was  measured  for  a  single  sine  wave  of  frequency  u>o  =  1.0  rad  I  sec  at 
tow  speeds  from  0.0  mj sec  to  5.0  m/ sec  at  a  tow  angle  of  45°.  The  resulting  plot  of  coe  vs.  U 
is  shown  in  Figure  2-16.  The  plot  of  u>e  vs.  U  is  expected  to  be  linear  with  a  slope  of  —  c°s^ 
(see  Equation  2.85),  which  for  =  1.0  rad/ sec  and  0  =  45°  has  a  value  of  -0.0721  m_1.  The 
measured  slope  was  -0.0767  m-1,  which  shows  the  proper  dependence  of  encounter  frequency 
on  speed  has  been  captured  by  the  model. 

2.5.3  Spatial  Spectral  Analysis 

By  sampling  the  wave  record  spatially  for  sea  surface  elevation,  one  can  determine  the  spatial 
power  spectral  density  curve  of  the  modeled  wave  field.  By  comparing  the  resultant  spectral 
density  curve  to  the  desired  spectral  density  curve,  a  verification  of  the  model’s  wave  field 
generation  method  can  be  conducted.  Figure  2-17  shows  the  spectral  analysis  of  spatial  wave 
data  for  a  wave  field  constructed  with  only  one  propagation  direction.  The  resultant  energy 
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Spatial  Power  Spectral  Density  Analysis 
3  Propagation  Directions 


Figure  2-18:  Spectral  Analysis  of  Spatial  Wave  Data  -  Three  Propagation  Directions 

density  spectrum  closely  resembles  the  desired  input  spectrum  and  there  is  good  agreement 
amongst  the  values  representing  the  area  under  the  input  spectrum,  the  area  under  the  output 
spectrum,  and  the  mean  square  wave  elevation  of  the  modeled  wave  field. 

An  interesting  phenomena  occurs  when  the  number  of  propagation  directions  is  increased 
to  three.  As  shown  in  Figure  2-18,  there  is  an  apparent  increase  in  energy  density  at  the  lower 
spatial  frequencies  while  the  opposite  occurs  at  the  higher  spatial  frequencies.  There  is  still 
good  agreement  amongst  the  values  representing  the  area  under  the  input  spectrum,  the  area 
under  the  output  spectrum,  and  the  mean  square  wave  elevation  of  the  modeled  wave  field, 
but  there  is  an  obvious  shift  in  energy  to  lower  frequencies.  This  shift  in  energy  to  lower 
frequencies  is  attributed  to  the  wave  summation  method  listed  in  Equation  2.73  and  shown  in 
Figure  2-15.  As  the  individual  wave  is  added,  it  is  being  ” stretched”  in  length  due  to  the  angle 
between  the  antenna  and  the  wave  itself.  For  the  extreme  case  where  this  angle  is  90°,  only  one 
value  of  wave  elevation  from  the  individual  wave  would  be  applied  to  the  entire  length  of  the 
antenna.  In  this  case,  a  single  point  value  is  stretched  along  the  entire  length  of  the  antenna. 
In  the  opposite  extreme  case  where  the  angle  is  0°,  there  would  be  a  point  to  point  match  in 
the  individual  wave  and  the  antenna  with  no  stretching,  as  shown  in  Figure  2-17.  The  same 
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Spatial  Power  Spectral  Density  Analysis 
5  Propagation  Directions 


Figure  2-19:  Spectral  Analysis  of  Spatial  Wave  Data  -  Five  Propagation  Directions 

concept  applies  to  angles  in  between  the  extreme  cases  as  a  small  portion  of  the  wave  field  is 
stretched  to  fit  the  antenna. 

The  stretching  causes  a  decrease  in  spatial  frequency  as  the  apparent  wavelength  of  the 
individual  wave  is  increased,  as  shown  in  Figure  2-20.  Since  the  individual  spatial  frequency 
decreases,  the  total  energy  of  the  wave  field  along  the  antenna  shifts  to  lower  frequencies.  It 
should  also  be  pointed  out  that  the  shifts  in  energy  magnitude  occur  at  values  that  correspond 
to  the  cosine  of  the  angle  between  the  individual  waves  and  the  antenna.  In  the  case  of  three 
propagation  directions,  shown  in  Figure  2-18,  the  waves  are  propagating  at  -60°,  0°,  and  +60°. 
The  shift  in  energy  magnitude  occurs  at  0.5  m_1,  which  is  the  same  value  as  the  width  of  the 
input  spectrum  multiplied  by  cos  (60°) .  In  the  case  of  five  propagation  directions,  shown  in 
Figure  2-19,  the  waves  are  propagating  at  -72°,  -36°,  0°,  +36°,  and  +72°.  The  shifts  in  energy 
magnitude  occur  at  0.3  m-1  and  0.8  m~l,  which  are  the  same  values  as  the  width  of  the  input 
spectrum  multiplied  by  cos  (36°)  and  cos  (72°),  respectively. 
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Figure  2-20:  Spatial  Frequency  Stretching 


2.5.4  Temporal  Spectral  Analysis 

By  sampling  the  wave  record  temporally  for  sea  surface  elevation,  one  can  determine  the  tem¬ 
poral  power  spectral  density  curve  of  the  modeled  wave  field.  By  comparing  the  resultant 
spectral  density  curve  to  the  desired  spectral  density  curve,  a  verification  of  the  model’s  wave 
field  generation  method  can  be  conducted.  Figure  2-21  shows  the  temporal  spectral  analysis 
of  a  model  generated  wave  field.  The  resultant  energy  density  spectrum  closely  resembles  the 
desired  input  spectrum  and  there  is  good  agreement  amongst  the  values  representing  the  area 
under  the  input  spectrum,  the  area  under  the  output  spectrum,  and  the  mean  square  wave 
elevation  of  the  modeled  wave  field. 


Temporal  Power  Spectral  Density  Analysis 


Figure  2-21:  Spectral  Analysis  of  Temporal  Wave  Data 
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Chapter  3 


Antenna  Model  Development 

3.1  Coordinate  System  Description 

The  coordinate  system  adopted  for  the  model  is  a  three-dimensional  system  whose  origin  trans¬ 
lates  in  the  antenna  tow  direction  at  the  antenna  tow  speed  U,  which  is  defined  to  be  positive. 
The  origin  of  the  antenna  is  defined  to  be  the  trailing  point  of  the  antenna.  The  waves  travel 
in  the  positive  x  direction,  all  at  varying  angles,  but  always  from  left  to  right  in  the  positive  x 
direction.  Figure  3-1  depicts  the  coordinate  system  in  the  x-z  plane  as  well  as  the  basic  antenna 
configuration.  The  antenna  tow  angle,  9,  is  the  angle  measured  from  the  antenna  to  the  prin¬ 
cipal  axis  of  wave  propagation.  It  is  defined  to  be  0°  when  pointing  in  the  positive  x  direction 
and  180°  when  pointing  in  the  negative  x  direction.  The  model  allows  the  user  to  select  any 
tow  angle  between  and  including  0°  and  180°.  Due  to  port  and  starboard  symmetry  of  the 
anntena,  the  behavior  of  the  antenna  for  the  angles  180°  to  360°  is  identical  to  the  behavior 
of  the  antenna  for  the  angles  0°  to  180°.  Figure  3-2  depicts  the  coordinate  system  in  the  x-y 
plane  as  well  as  the  antenna  tow  direction. 

The  equations  listed  in  Table  2.1  apply  to  a  stationary  coordinate  system.  The  transforma¬ 
tion  to  a  moving  coordinate  system  employs  the  following  relationship 


x'  =  x  +  Ut 


(3.1) 


where 
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Direction  of 
Wave  Propagation 


Tow  Point 


Figure  3-1:  Coordinate  System:  X-Z  Plane 


Y 


Head 

Seas 
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Seas 


Figure  3-2:  Coordinate  System:  X-Y  Plane 


x'  is  the  earth-fixed,  stationary,  coordinate  system 

x  is  the  moving  coordinate  system  whose  origin  translates  with  velocity  U  ,  as  shown  in 
Figures  3-1  and  3-2. 


3.2  Governing  Vertical  Equation  of  Motion 

The  vertical  forces  acting  on  the  antenna  are  given  in  the  right  hand  side  of  Equation  3.2 

m<^t =  T('x’  ^^+Cd^w~  + b(x’  + Cma  pwater  717,2  + (3-2) 


where 

m  is  the  mass  per  unit  length  of  the  antenna,  assumed  uniform 
v  is  antenna  vertical  velocity,  which  is  equal  to  ^ 

T(x,t)  is  the  tension  along  the  length  of  the  antenna,  arising  from  skin  frictional  drag 
between  the  submerged  portion  of  the  antenna  and  the  water 

q  is  the  z  coordinate  of  the  antenna  central  axis,  as  shown  in  Figure  3-3 
Cd  is  the  normal  drag  coefficient  of  a  cylindrical  antenna  section  as  it  moves  in  the  z 
direction  through  the  water 

w  is  the  fluid  vertical  velocity 

b(x,t)  is  the  net  buoyancy  force  acting  on  the  antenna,  which  is  the  difference  between  the 
downward  force  due  to  its  weight  and  the  upward  force  due  to  the  displaced  water. 

Cma  is  the  coefficient  of  added  mass  for  a  cylindrical  section  in  a  flow  field  normal  to  its 
long  axis 

Pwater  IS  the  density  of  water  (for  seawater  of  15°  C,  pwater  —  1025.9  kg/rv? ) 
r  is  the  antenna  radius 

El  is  the  product  of  the  antenna  cross  sectional  modulus  of  elasticity  and  area  moment  of 
inertia,  thus  a  measure  of  lateral  stiffness. 

The  first  term  on  the  right  hand  side  of  Equation  3.2  represents  a  restoring  force  charac¬ 
teristic  of  a  string  displaced  from  its  equilibrium  condition.  The  tension  in  the  antenna  is 
calculated  from  the  time  varying  wetted  surface  area  along  the  length  of  the  antenna.  A  tan¬ 
gential  drag  coefficient  of  Ct  =  0.0035  was  used  based  upon  Reynold’s  numbers  which  ranged 
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z 


Figure  3-3:  Antenna  and  Sea  Surface  Elevation  Variables 

from  32,940  to  592,900  for  speed/diameter  combinations  of  3  knots/1  inch  to  9  knots/6  inches, 
respectively.  [1]  The  second  term  represents  the  normal  drag  force  a  cylindrical  element  would 
experience  as  it  is  moved  through  the  fluid  with  its  long  axis  perpendicular  to  the  direction  of 
motion.  A  value  of  Cp  =  1.0  was  used,  which  assumes  an  infinite  fluid  domain.  The  fourth 
term  represents  the  added  mass  effect  experienced  by  a  cylindrical  object  as  the  fluid  around  it 
is  accelerating  normal  to  the  long  axis  of  the  cylinder.  A  value  of  Cm  A  =  1.0  was  used,  which 
is  consistent  with  standard  fluid  dynamic  texts  assuming  an  infinite  fluid  domain.  [6]  The  final 
term  represents  the  restoring  force  due  to  the  stiffness  of  the  material,  and  its  distribution  with 
respect  to  the  neutral  axis.  Initial  design  options  have  the  lateral  stiffness  ranging  from  300  to 
1400  Pa-rn2. 

It  is  recognized  that  some  of  these  forces  have  been  approximated,  namely  in  that  the  values 
quoted  above  are  for  Cp  and  Cm  A  in  an  infinite  fluid  domain.  The  antenna  operates  in  close 
proximity  to  the  free  surface,  and  occasionally  portions  rise  completely  above  the  surface  of 
the  water.  A  rigorous  treatment  of  the  behavior  of  Cp  and  Cm  A  near  the  free  surface  was  not 
attempted.  Rather,  a  straightforward  linear  approximation  was  employed  which  varied  the 
coefficients  Cp  and  Cm  A  from  their  infinite  fluid  domain  values  when  the  antenna  axis  was  at 
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Conductor: 
Isometric  View 


Antenna  Element: 
Section  View 


Figure  3-4:  Antenna  Element  Description 

a  depth  of  r,  just  submerged,  to  a  value  of  zero  when  the  antenna  axis  was  at  a  height  of  r, 
just  breaking  out  of  the  water.  This  dependence  of  Cd  and  Cm  A  on  relative  antenna  depth,  y, 
is  shown  in  Equations  3.3  and  3.4,  which  only  apply  to  the  partially  submerged  case. 


Cd{v)  =  CDlnfinite{-^  +  i) 

(3.3) 

y  i 

Cma(v)  =  CMAInfinite(-7^  +  y) 

(3.4) 

3.3 

Antenna  Modeling  Methods 

3.3.1 

Discretization 

The  antenna,  of  length  L,  is  subdivided  first  into  lengths  of  the  physical  radio  antenna  elements. 
These  antenna  elements,  designed  at  Lincoln  Laboratory,  are  called  “jelly  roll”  antennas  because 
they  are  a  conducting  sheet  rolled  up  about  1^  times.  A  graphical  description  of  the  antenna 
element  is  shown  in  Figure  3-4.  The  element  length  is  designated  Leiem.  This  first  level  of 
subdivision  is  not  only  logical,  but  necessary  for  the  exposure  statistics  calculated  by  the  model. 
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Figure  3-5:  Antenna  and  Antenna  Element  Description 

As  it  turns  out,  this  type  of  antenna  element  performs  well  as  long  as  it  has  some  small  amount, 
on  the  order  of  1  mm,  of  its  outer  housing  diameter  exposed  along  its  entire  length.  But  if 
the  element  becomes  covered  by  even  a  small  amount  of  water  at  any  point  along  its  length, 
antenna  performance  drops  dramatically.  Thus  it  is  necessary  to  monitor  numerous  points 
within  each  element  to  determine  whether  a  given  element  is  exposed  or  not.  This  leads  to 
the  final  level  of  antenna  subdivision:  each  element  is  meshed  into  some  integer  number  M  of 
sections  with  length  As.  The  subdivision  process  is  shown  in  Figure  3-5. 

An  interesting  feature  of  the  models  which  calculate  the  sea  elevation  and  fluid  velocities 
by  Fast  Fourier  Transforms  (FFTs)  is  that  the  frequency  meshing  (Aw)  of  the  spectral  energy 
density  function  S(u),  is  now  related  to  the  spatial  meshing  (Ax)  of  the  sea,  and  thus  the 
spatial  meshing  of  the  antenna  (As).  This  means  that  there  is  an  extra  constraint  to  consider 
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in  selecting  the  finest  antenna  mesh  size  As;  or  conversely,  given  a  mesh  size  As,  constraints 
are  imposed  on  the  Au;  meshing  of  the  spectral  energy  density  function  S(u>).  Some  freedom  is 
regained  if  one  decouples  the  length  of  the  antenna  and  the  length  of  the  sea  simulated.  This 
decoupling  is  implemented  in  the  model,  so  that  for  a  typical  run  the  antenna  may  be  25  m 
long  with  a  mesh  size  of  0.0625  m,  while  the  length  of  sea  modeled  is  512  m  long  with  a  mesh 
size  of  0.0625  m.  The  key  point  is  that  the  mesh  length  between  the  antenna  (As)  and  the  sea 
(Ax)  must  match. 

3.3.2  Integration  Method 

The  time  marching  integration  method  implemented  is  an  explicit  forward  Euler  scheme.  The 
following  discussion  of  of  the  explicit  forward  Euler  scheme  comes  from  Ulrich.  [9]  Such  schemes 
are  characterized  by  using  information  at  time  t  to  calculate  functions  at  time  t  +  At.  A  simple 
example  of  how  the  forward  Euler  method  is  used  is  given  below.  Consider  the  following  sample 
numerical  problem. 

Solve 


du 

dt 

u(0) 


=  u 


=  1 


The  time  derivative  is  approximated  using  the  Euler  method  as 

Un+ 1  -  Un 


du 

~dt 


t=tn 


which  for  constant  time  step  A t  becomes 


du 

dt 


t—tn 


tn+ 1  -  tn 


jjn+l  _  jjn 

At 


having  adopted  the  notation  of  Celia  and  Gray  [2] ,  where 


(3.5) 

(3.6) 


(3.7) 


(3.8) 


Un  =  u{tn ) 


(3.9a) 
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and  tn  is  the  nth  time  step.  Then  the  solution  proceeds  step  by  step,  beginning  with  the  initial 
value  of  u  given  in  Equation  3.6 


du 

dt 


=  -u(O)2  =  1 


(3.10) 


t= o 


then  solving  Equation  3.8  for  Un+1  we  obtain 


U1  =  U°  +  At  ^ 
dt 


(3.11) 


t= o 


With  the  value  of  u  now  known  at  time  t  (f/1),  we  can  calculate 


du 

dt 


U2 


t- 1 


=  ([/■) 


(3.12) 


and 


(3.13) 

In  like  manner,  we  move  forward  in  time  to  calculate  any  desired  Un.  Numerical  experi¬ 
mentation  was  used  to  determine  the  appropriate  time  step  size  which  produced  acceptable 
results. 

The  spatial  derivatives  (^)  used  in  the  model  are  classified  as  centered  finite-divided- 
difference  formulas.  [3]  Thus  the  antenna  algorithm  can  be  classified  as  a  forward  time,  centered 
space  (FTCS)  method. 

3.3.3  Boundary  Conditions 

The  boundary  conditions  imposed  on  the  antenna  endpoints  are  given  in  Equations  3.14  and 
3.15  as 


U2  =  U1  +  At  ^ 
dt 


q(  o) 

=  Mo) 

(3.14) 

q(L) 

=  h(L) 

(3.15) 
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That  is,  the  antenna  follows  the  sea  surface  at  its  endpoints.  This  is  envisioned  to  be  accom¬ 
plished  at  the  trailing  end  by  a  drogue  having  positive  buoyancy,  and  at  the  tow  point  by  one 
of  several  options.  These  options  include 


•  a  dynamic  lifting  body 

•  a  surface  following  buoy 

•  a  leader  of  positive  buoyancy 

•  some  combination  of  the  above  features. 

3.3.4  Initial  Conditions 

The  antenna  simulation  begins  at  t  =  0  with 

q(x)  —  h(x) 

dq{x)  _  dh(x ) 
dt  dt 


(3.16) 

(3.17) 


that  is,  the  antenna  is  placed  at  the  sea  surface,  and  given  an  initial  vertical  velocity  matching 
that  of  the  sea  surface. 
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Chapter  4 


Modeling  Results 

4.1  Exposure  Statistics 

4.1.1  Definition  of  Exposure 

The  primary  goal  of  the  modeling  effort  was  to  predict  exposure  statistics  as  various  design 
and  operational  parameters  were  varied;  therefore,  it  is  necessary  to  define  exposure  before 
proceeding  on  to  the  exposure  statistics  results.  The  overall  exposure  of  the  antenna  is  driven 
by  the  exposure  of  individual  antenna  elements,  so  the  determination  of  an  exposed,  or  receiving, 
antenna  element  must  be  addressed.  The  experimentally  observed  performance  characteristics 
of  the  “jelly  roll”  antenna  elements  led  to  the  following  definition  of  exposure.  It  was  found  that 
for  a  level  element  in  calm  water,  the  antenna  gain  is  stable  right  up  to  the  point  at  which  the 
outer  skin  becomes  completely  submerged.  It  was  also  found  that  if  any  portion  of  the  outer 
skin  was  completely  covered  by  water  then  the  gain  of  the  element  suffered  drastically.  Thus  an 
exposed  element  was  defined  to  be  one  which  was  exposed  by  some  small  threshold  height,  set 
to  four  percent  of  Diameter  for  all  simulations,  at  every  point  along  the  length  of  the  element. 
In  a  discrete  computer  model  some  approximation  is  necessary  with  regard  to  every  point  along 
the  length  of  the  element.  Thus  each  antenna  element  was  divided  into  smaller  segments  for 
analysis.  Typical  designs  considered  had  an  element  length  of  0.5  m  and  were  further  divided, 
or  oversampled,  at  M  —  8. 
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4.1.2  Statistics  Calculated  by  the  Model 

The  computer  model  was  programmed  to  track  three  key  statistics: 

•  average  exposure  -  the  average  fraction  of  antenna  elements  exposed  throughout  the  entire 
duration  of  the  simulation 

•  threshold  met  -  the  average  fraction  of  the  time  that  some  specified  number  of  antenna 
elements,  or  more,  were  exposed 

•  average  counterdetection  length  -  the  average  length  of  the  antenna  where  there  is  no 
water  touching  the  antenna. 

In  addition,  a  record  is  made  at  each  sample  time  of  the  condition  of  every  antenna  element, 
such  that  when  the  simulation  is  finished,  an  exposure  matrix  is  available  for  further  post¬ 
processing.  Examples  of  information  which  could  be  extracted  from  the  exposure  matrix 
include 

•  average  drop-out  time  of  a  specific  element 

•  average  spatial  length  of  wave  wash-overs 

•  histograms  and  standard  deviations  of  the  number  of  exposed  elements 

•  exposure  statistics  for  an  array  not  fully  populated  (i.e.,  one  could  examine  only  the 
columns  which  had  antenna  elements  installed)  . 

Five  output  files  are  generated  for  each  run  of  the  computer  model.  These  output  files  are 
described  in  detail  in  Appendix  A.  The  input  parameters  which  are  specified  for  each  run  of 
the  computer  model  are  also  described  in  detail  in  Appendix  A. 

In  the  exposure  statistics  results  given  below,  each  series  of  simulations  is  presented  in  the 
following  manner: 

•  a  listing  of  the  antenna  and  operational  parameters  input  table  used  to  generate  the  series 

•  a  plot  of  the  average  fraction  of  elements  exposed 
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•  a  brief  discussion  of  the  exposure  plot. 

It  is  important  to  realize  that  these  plots  represent  averages  in  time  and  space.  This  is 
necessary  in  order  to  display  the  results  as  functions  of  key  antenna  design  and  operational 
parameters,  which  is  most  useful  from  a  designer’s  viewpoint.  From  the  output  files,  it  is  also 
possible  to  analyze  how  the  statistics  vary  as  time  progresses,  or  how  the  statistics  vary  along 
the  array  from  the  tow  point  to  the  trailing  end.  Although  there  are  a  large  number  of  possible 
statistics  that  could  be  analyzed  for  an  even  larger  number  of  design  and  operational  param¬ 
eter  combinations,  only  one  statistic  is  presented  here  for  seven  selected  design  or  operational 
parameters. 
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Parameter  Description 


Duration 

Nelern 

Lant 

Lelem 

Diameter 

El 

Pant 

M 

U 

Tow  Angle 
Hs 

Development 
Spectral  Model 
Spreading  Function 
Nfreq 

NN 


60  sec 
73 

36.5  m 
0.5  m 
3.0  inches 
100  N-m2 
408  kg/m3 
8 

5.0  kts 

0°  ,45°  ,90°  ,135° ,  1 80' 
3.0  m 
1.0 

Bretschneider 

cos2 


213 

1,5,9,13,17,21 


Table  4.1:  Series  1  Input  Parameters 


4.2  Results 

4.2.1  Series  1:  Statistics  vs.  Number  of  Propagation  Directions 
Input 

The  input  parameters  for  Series  1:  Statistics  versus  Number  of  Propagation  Directions,  are 
given  in  Table  4.1.  This  is  a  measure  of  the  fineness  of  the  angular  sea  description.  The 
number  of  propagation  directions  used  were  1,  5,  9,  13,  17,  and  21.  The  tow  angles  used  were 
0,  45,  90,  135,  and  180  degrees. 


Exposure  \s  Number  of  Propagation  Directions 


Figure  4-1:  Average  Exposure  vs.  Number  of  Propagation  Directions 


Results 

Figure  4-1  shows  a  very  weak  dependence  of  the  average  fraction  of  elements  exposed  on  the 
number  of  propagation  directions  used  by  the  model  to  simulate  the  multidirectional  wave  field. 
However,  the  plot  shows  a  strong  dependence  of  average  exposure  on  tow  angle.  The  insensi¬ 
tivity  to  the  number  of  propagation  directions  in  the  sea  description  not  only  demonstrates  the 
validity  of  the  approach  used,  but  also  shows  that  computing  time  can  be  minimized  by  using 
just  a  modest  number  (3  to  9)  of  directions  in  the  sea  description. 
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Parameter 

Description 

Duration 

60  sec 

Nelem 

73 

Lant 

36.5  m 

Lelem 

0.5  m 

Diameter 

3.0  inches 

EI 

100  N-m 2 

Pant 

408  kg/m 3 

M 

8 

U 

5.0  kts 

Tow  Angle 

0°, 30°, 60°  ,90°, 120°, 150°, 180° 

Hs 

3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function  cos2 

Nfreq 

213 

NN 

3 

Table  4.2:  Series  2  Input  Parameters 


4.2.2  Series  2:  Statistics  vs.  Tow  Angle 


Input 

The  input  parameters  for  Series  2:  Statistics  versus  Tow  Angle,  are  given  in  Table  4.2.  The 
tow  angles  used  were  0,  30,  60,  90,  120,  150,  and  180  degrees. 
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Exposure  vs  Tow  Angle 


Figure  4-2:  Average  Exposure  vs.  Tow  Angle 

Results 

Figure  4-2  shows  a  rather  weak  dependence  of  the  average  fraction  of  elements  exposed  on  the 
tow  angle  from  0°-90°  (following  seas);  however,  there  is  a  very  strong  dependence  on  tow  angle 
from  90°-180°  (head  seas). 
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Parameter 

Description 

Duration 

60  sec 

N elem 

73 

Lant 

36.5  m 

Lelem 

0.5  m 

Diameter 

3.0  inches 

EI 

100  N-m2 

Pant 

408  kg/m 3 

M 

8 

U 

5.0  kts 

Tow  Angle 

0°, 45°, 90°, 135°, 180° 

Hs 

0.5, 1.0, 1.5, 2.0, 2.5, 3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function 

cos2 

Nfreq 

213 

NN 

3 

Table  4.3:  Series  3  Input  Parameters 


4.2.3  Series  3:  Statistics  vs.  Sea  Severity 
Input 

The  input  parameters  for  Series  3:  Statistics  versus  Sea  Severity,  are  given  in  Table  4.3.  The 
significant  wave  heights  used  were  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  meters.  The  tow  angles  used 
were  0,  45,  90,  135,  and  180  degrees. 


Exposure  vs  Significant  Wa\«  Height 


Figure  4-3:  Average  Exposure  vs.  Sea  Severity 


Results 

Figure  4-3  shows  a  rather  weak  dependence  of  the  average  fraction  of  elements  exposed  on  the 
significant  wave  height  at  all  tow  angles. 
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Parameter 

Description 

Duration 

60  sec 

Nelem 

73 

Lant 

36.5  m 

Lelem 

0.5  m 

Diameter 

3.0  inches 

EI 

200  N-m2 

Pant 

408  kg/m 3 

M 

8 

U 

3, 5, 7, 9  kts 

Tow  Angle 

0°, 45°  ,90°, 135°, 180° 

Hs 

3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function 

cos2 

Nfreq 

213 

NN 

3 

Table  4.4:  Series  4  Input  Parameters 


4.2.4  Series  4:  Statistics  vs.  Tow  Speed 
Input 

The  input  parameters  for  Series  4:  Statistics  versus  Tow  Speed,  are  given  in  Table  4.4.  The 
tow  speeds  used  were  3,  5,  7,  and  9  knots.  The  tow  angles  used  were  0,  45,  90,  135,  and  180 
degrees. 
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Exposure  vs  Tow  Speed 


Figure  4-4:  Average  Exposure  vs.  Tow  Speed 


Results 

Figure  4-4  demonstrates  the  effect  that  tow  angle  has  on  other  parameters.  In  this  instance, 
increasing  tow  speed  for  tow  angles  greater  than  90°  (head  seas)  causes  a  decrease  in  exposure 
fraction.  While,  for  tow  angles  less  than  90°  (following  seas),  increasing  tow  speed  causes  an 
increase  in  exposure  fraction. 
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Parameter 

Description 

Duration 

60  sec 

Nelem 

10,30,50,120,200 

Lant 

5,15,25,60,100  m 

Lelem 

0.5  m 

Diameter 

3.0  inches 

EI 

200  N-m 2 

Pant 

408  kg/m 3 

M 

8 

U 

5.0  kts 

Tow  Angle 

0°, 45°  ,90°, 135°, 180° 

Hs 

3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function  cos2 

Nfreq 

213 

NN 

3 

Table  4.5:  Series  5  Input  Parameters 


4.2.5  Series  5:  Statistics  vs.  Antenna  Length 
Input 

The  input  parameters  for  Series  5:  Statistics  versus  Antenna  Length,  are  given  in  Table  4.5. 
The  antenna  lengths  used  were  5,  15,  25,  60,  and  100  meters.  The  tow  angles  used  were  0,  45, 
90,  135,  and  180  degrees. 
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Exposure  vs  Antenna  Length 


Figure  4-5:  Average  Exposure  vs.  Antenna  Length 


Results 

Figure  4-5  shows  a  very  strong  dependence  of  the  average  fraction  of  elements  exposed  on  the 
antenna  length  for  lengths  less  than  25  meters.  For  lengths  greater  than  25  meters ,  there  is 
little  or  no  dependence  on  antenna  length. 
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Parameter 

Description 

Duration 

60  sec 

Nelem 

73 

L>ant 

36.5  m 

-Lelem 

0.5  m 

Diameter 

3. 0,6.0  inches 

El 

200  N-m 2 

Pant 

408  kg/m 3 

M 

8 

U 

5.0  kts 

Tow  Angle 

0°  ,45°  ,90°  ,135°  ,180° 

HS 

3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function 

0 

COS^ 

Nfreq 

213 

NN 

3 

Table  4.6:  Series  6  Input  Parameters 


4.2.6  Series  6:  Statistics  vs.  Antenna  Diameter 
Input 

The  input  parameters  for  Series  6:  Statistics  versus  Antenna  Diameter,  are  given  in  Table  4.6. 
The  antenna  diameters  used  were  3  and  6  inches.  The  tow  angles  used  were  0,  45,  90,  135, 
and  180  degrees. 
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Exposure  vs  Antenna  Diameter 


Figure  4-6:  Average  Exposure  vs.  Antenna  Diameter 


Results 

Figure  4-6  also  demonstrates  the  influence  that  tow  angle  has  on  amplifying  the  effects  of  design 
parameter  changes.  In  all  cases  shown,  the  larger  diameter  antenna  has  an  increased  exposure 
fraction,  but  the  change  is  most  significant  for  head  seas  cases. 
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Parameter 

Description 

Duration 

60  sec 

Nelem 

73 

Lant 

36.5  m 

Lelem 

0.5  m 

Diameter 

3.0  inches 

EI 

200  N-m2 

Pant 

200,300,400,500,600  kg/m3 

M 

8 

U 

5.0  kts 

Tow  Angle 

0°,  45°,  90°,  135°,  180° 

Hs 

3.0  m 

Development 

1.0 

Spectral  Model 

Bretschneider 

Spreading  Function 

cos2 

Nfreq 

213 

NN 

3 

Table  4.7:  Series  7  Input  Parameters 


4.2.7  Series  7:  Statistics  vs.  Antenna  Density 
Input 

The  input  parameters  for  Series  7:  Statistics  versus  Antenna  Density,  are  given  in  Table  4.7. 
The  antenna  densities  used  were  200,  300,  400,  500,  and  600  kg/m3.  The  tow  angles  used  were 
0,  45,  90,  135,  and  180  degrees. 
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Exposure  vs  Antenna  Density 


Figure  4-7:  Average  Exposure  vs.  Antenna  Density 


Results 

Figure  4-7  again  demonstrates  the  amplifying  effect  that  tow  angle  has  on  design  parameter 
changes.  In  all  cases  shown,  as  you  increase  the  antenna  density  there  is  a  resulting  decrease  in 
exposure  fraction;  however,  this  decrease  is  amplified  for  tow  angles  in  the  head  seas  direction. 
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Chapter  5 

Future  Efforts 


The  combined  Navy  and  DARPA  Multi-element  Buoyant  Cable  Antenna  effort  was  funded 
in  FY  00  for  3  years.  The  demonstration  phase  culminated  with  a  successful  testing  effort  in 
Hawaii  in  the  Summer  of  1999.  Within  FY  00,  Lincoln  Laboratory.is  scheduled  to  complete  the 
development  and  design  of  an  antenna  which  can  stand  up  to  the  rigors  of  an  at  sea  deployment. 
FY  01  will  see  the  initial  production  of  the  deployable  version  with  the  goal  of  at  sea  submarine 
testing  during  FY  02. 

Work  at  MIT  will  continue  to  support  the  buoyant  cable  antenna  effort.  Another  graduate 
student,  under  the  supervision  of  Professor  Jerome  Milgram,  will  continue  and  extend  the 
scope  of  the  simulation  effort.  Areas  which  will  be  investigated  to  improve  the  capability  of 
the  computer  model  include 

•  more  detailed  modeling  of  the  hydrodynamics  at  the  free  surface 

•  higher  order  integration  schemes 

•  addition  of  digital  signal  processors  to  improve  computation  speed 

•  validation  against  tow  tank  data 

•  validation  against  at-sea  test  data 

•  effects  of  roll  for  revised  antenna  elements  whose  performance  depends  on  roll  angle. 
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Appendix  A 


User’s  Guide  to  Proteus 

A.l  Introduction 

A.  1.1  Contents  of  the  User’s  Guide 

The  User’s  Guide  is  broken  down  into  three  sections: 

•  Introduction 

•  Input  File 

•  Output  Files. 

Section  A.l  describes  the  contents  of  the  user’s  guide  and  contains  the  procedures  necessary 
to  obtain  the  results  discussed  in  Section  A. 3.  Section  A. 2  describes  the  Proteus  input  file  line 
by  line  and  includes  program  features  and  limitations  associated  with  each  parameter.  Section 
A. 3  briefly  describes  the  contents  of  each  of  the  five  output  files  and  the  recommended  methods 
of  viewing  or  extracting  the  information. 

A.  1.2  Analysis  Procedure 
Obtaining  the  Snapshot  Output 

1.  Ensure  that  the  input  file,  proteus_input.txt,  is  located  in  the  same  directory  as  the 
FORTRAN  executable  file,  proteus.  exe. 
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2.  Modify  the  input  file  to  reflect  the  conditions  that  the  user  wishes  to  analyze,  using  Section 
A.2  as  guidance  for  features  and  limitations.  Ensure  that  the  matlab_  output  variable  is 
set  to  .TRUE. 

3.  Execute  proteus.  exe. 

4.  The  snapshot  output  file  (run_id/ _ mat’)  will  be  located  in  the  same  directory  as  the 
executable  file. 

5.  Ensure  that  the  output  file  {run_id/ _ mat’)  will  be  located  in  the  same  directory  as 
the  Matlabw  function  file,  proteus _  mat. 

(fD 

6.  At  the  Matlat)  command  line  type:  clear;  proteus _  mat(’run_  id/ / _  mat  ’). 

Obtaining  the  Parameters  and  Results  Summary 

1.  Ensure  that  the  input  file,  proteus_input.txt,  is  located  in  the  same  directory  as  the 
FORTRAN  executable  file,  proteus.exe. 

2.  Modify  the  input  file  to  reflect  the  conditions  that  the  user  wishes  to  analyze,  using  Section 
A.2  as  guidance  for  features  and  limitations.  Ensure  that  the  run_  info  variable  is  set  to 
.TRUE. 

3.  Execute  proteus.exe. 

4.  The  parameters  and  results  summary  output  file  {run_id/ _ info.txt’)  will  be  located  in 
the  same  directory  as  the  executable  file. 

5.  The  output  file  can  be  viewed  with  any  standard  text  editor. 

Obtaining  the  Exposure  Data  Matrix 

1.  Ensure  that  the  input  file,  proteus_input.txt,  is  located  in  the  same  directory  as  the 
FORTRAN  executable  file,  proteus.exe. 

2.  Modify  the  input  file  to  reflect  the  conditions  that  the  user  wishes  to  analyze,  using  Section 
A.2  as  guidance  for  features  and  limitations.  Ensure  that  the  exp _  data  variable  is  set  to 
.TRUE. 
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3.  Execute  proteus.exe. 


4.  The  exposure  data  matrix  output  file  (run_id/ /' _e dm’)  will  be  located  in  the  same 
directory  as  the  executable  file. 


5. 


(K) 

The  output  file  can  be  viewed  with  a  user  supplied  Matlao^  file  or  with  a  spreadsheet 
application. 


Obtaining  Wave  Data  Temporal  Spectrum  Analysis 

1.  Ensure  that  the  input  file,  proteus_input.txt,  is  located  in  the  same  directory  as  the 
FORTRAN  executable  file,  proteus.exe. 

2.  Modify  the  input  file  to  reflect  the  conditions  that  the  user  wishes  to  analyze,  using  Section 
A.2  as  guidance  for  features  and  limitations.  Ensure  that  the  wave_  temporal  variable  is 
set  to  .TRUE. 

3.  Execute  proteus.exe. 

4.  The  temporal  wave  data  output  file  {run_id/  f  _ mat2’)  will  be  located  in  the  same  di¬ 
rectory  as  the  executable  file. 

5.  Ensure  that  the  output  file  (run_id/ /’ _ mat2’)  will  be  located  in  the  same  directory  as 
the  executable  file,  proteus _  mat2.  exe. 

6.  Execute  proteus_mat2.exe.  When  the  program  prompts  the  user  for  filename,  enter 
run_  id/ '/ _  mat2. 

7.  Ensure  that  the  output  file  (run_id/ /, _ mat2_proc’)  will  be  located  in  the  same  directory 

(?) 

as  the  MaMaf)  function  file,  proteus_  mat2. 

(R) 

8.  At  the  Matlabw  command  line  type:  clear;  proteus_mat2(’run_id// _mat2_proc’). 

Obtaining  Wave  Data  Spatial  Spectrum  Analysis 

1.  Ensure  that  the  input  file,  proteus_input.txt,  is  located  in  the  same  directory  as  the 
FORTRAN  executable  file,  proteus.exe. 
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2.  Modify  the  input  file  to  reflect  the  conditions  that  the  user  wishes  to  analyze,  using  Section 
A.2  as  guidance  for  features  and  limitations.  Ensure  that  the  wave_  spatial  variable  is 
set  to  .TRUE. 


3.  Execute  proteus,  exe . 

4.  The  spatial  wave  data  output  file  ( run__id / will  be  located  in  the  same  directory 
as  the  executable  file. 


5.  Ensure  that  the  output  file  (run_id/ / _ mat35)  will  be  located  in  the  same  directory  as 
the  executable  file,  proteus _  mats.  exe. 


6.  Execute  proteus_mat3.exe.  When  the  program  prompts  the  user  for  filename,  enter 
run_  id/ /__  mats . 


7. 


Ensure  that  the  output  file  (rwi_Jd//’_mat3_proc’)  will  be  located  in  the  same  directory 
(K) 

as  the  Matlat)  function  file,  proteus _  mat3. 


8. 


At  the  Matlab 


command  line  type: 


clear;  proteus __  mat3(’run_  id/ /_  mat3_  proc  ’). 


A. 2  Proteus  Input  File 

A. 2.1  Description  of  Features  and  Limitations 
Header 

The  input  file  for  Proteus  is  shown  in  Section  A.2.2.  The  first  7  lines  are  for  informational 
purposes  and  have  no  effect  on  the  program.  However,  note  that  all  variables  are  given  in  SI 
units  and  that  there  must  be  7  lines  prior  to  the  first  variable  defined. 

Run  Identifier  ( run_  id) 

The  run_  id  is  an  8  character  variable  used  to  name  the  5  output  files.  The  5  output  files  are 
named 

•  run_id/ /'  _  mat’ 
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•  run_  id//’’  _info.txt’ 


•  run_id/ /’  _e dm’ 

•  run_id/ _  mat2’ 

•  run_id/ /’  _  mat3’. 

(The  / /  represents  FORTRAN90  string  concatenation.)  Output  files  are  located  in  the  local 
directory  from  which  the  program  is  run.  A  sample  pattern  for  run_id  is  run_ id=’15May_01\ 

(S) 

Characters  which  Matlah*  would  interpret  as  operators  (/,  +,  *,  etc.)  should  be  avoided  in 

(P) 

run_id  if  one  of  the  Matlabw  output  options  is  selected. 

FFT  Power  of  2  ( nnf ) 

The  FFT  algorithm  is  designed  to  operate  on  arrays  which  have  an  integer  power  of  2  elements. 
nnf  sets  this  power  of  2,  making 


N  sea  =  2nnf 


(A.l) 


where 

N_sea  is  the  number  of  sea  frequency  components  and  spatial  locations. 

The  length  of  the  sea  surface  modeled  in  the  spatial  domain  is  generally  much  larger  than 
the  antenna  length  requires;  however,  the  speed  advantage  of  the  FFT  method  more  than 
outweighs  the  penalty  of  carrying  along  the  extra  sea  information. 

Number  of  Antenna  Elements  ( num_  of_  elem) 

This  is  the  number  of  physical  antenna  receiver  elements  in  the  array,  assuming  that 

•  the  array  is  fully  populated 

•  the  elements  have  zero  space  between  them. 

Antenna  Element  Length  (element  length) 

This  is  the  length  of  each  antenna  element. 
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Antenna  Oversample  Factor  (M) 

Each  antenna  element  is  broken  down  into  an  integer  number,  M,  of  sections.  These  sections 
correspond  to  the  antenna  spatial  mesh  size  As.  By  enforcing  the  condition  that  As  =  Ax,  the 
computer  model  constructs  the  spatial  frequency  meshing,  Ab,  such  that 


h  1 

2  As 

(A.2) 

L_sea  =  As  N_sea 

(A.3) 

Ab  =  -A— 

L  sea 

(A.4) 

By  invoking  the  dispersion  relation,  we  convert  from  the  b  domain  to  the  u  domain,  using 

lj  =  \/2ngb  (A.5) 

This  process  shows  how  the  spatial  meshing  relates  to  the  sea  energy  frequency  domain 
meshing.  As  one  can  see,  the  fact  that  N_sea  is  constrained  to  be  a  power  of  2  introduces 
some  limitations  on  the  mesh  size  As,  since  it  is  required  that 

element _length  =  M  As  (A. 6) 

Tow  Speed  (  U ) 

The  tow  speed  of  the  antenna  is  assumed  constant  with  respect  to  time.  Tow  speed  must  be 
greater  than  zero  (in  order  to  generate  tension).  When  9,  the  angle  between  an  individual 
propagation  direction  and  the  antenna  tow  angle,  is  less  than  90  degrees  (i.e.  following  seas), 
sea  frequency  components  above  omega_crit, 

omega  crit  =  —A  -  (A.7) 

~~  UcosO 

are  zeroed  out  to  prevent  antenna  folding.  If  omega _crit  is  less  than  1.5  rad/sec,  then  a 
warning  is  issued  that  excessive  frequency  components  are  being  zeroed  out,  and  program 
execution  stops. 
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Simulation  Time  ( duration ) 

The  total  length  of  time  simulated  by  the  program. 

Antenna  Time  Step  Increment  (  t_  incr) 

The  time  step  increment  used  in  solving  the  antenna’s  governing  equation.  A  typical  value 
for  an  antenna  with  no  stiffness  is  0.001  seconds.  When  stiffness  is  added,  the  time  increment 
should  be  reduced  to  0.0001  seconds.  In  all  cases,  if  the  program  stops  executing  due  to  a 
run-time  math  error,  the  time  increment  should  be  made  smaller  to  allow  proper  execution  of 
the  program. 

Snapshot  Output  Interval  ( pic_time ) 

This  variable  allows  the  user  to  control  the  time  interval  between  snapshots,  which  are  recorded 
to  file  run_id/ mat’.  These  snapshots  record  both  the  sea  and  antenna  elevation  data  along 
the  entire  length  of  the  antenna,  at  a  resolution  of  As. 

Time  of  First  Snapshot  ( t_first_pic ) 

The  time  to  begin  recording  snapshots  is  user  defined.  This  is  useful  to  ignore  transient  behavior 
at  the  start  of  the  simulation,  or  to  save  only  a  representative  sample  at  the  end  of  a  lengthy 
simulation. 

Statistics  Output  Interval  ( stat_time ) 

This  variable  allows  the  user  to  control  the  time  interval  between  exposure  data  matrix  records, 
which  are  recorded  to  file  run_id/ /’ _e dm’.  Selection  of  an  excessively  small  stat_  time  results 
in  subsequent  records  containing  essentially  the  same  information,  and  contributes  little  to  the 
statistics  of  the  sample. 

Snapshot  Output  Option (matlab_  output ) 

The  snapshot  output  file  (run_id//'_ mat’)  is  optional,  since  it  does  not  contain  any  unique 
statistical  information,  and  requires  a  relatively  large  amount  of  disk  space.  The  primary 
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benefit  of  examining  the  snapshot  output  file  is  to  learn  the  physical  behavior  of  the  antenna  in 
a  seaway.  In  addition,  the  output  is  useful  to  detect  when  the  time  step  size  is  too  small,  which 
is  characterized  by  rapid  unphysical  antenna  oscillations,  matlab  _  output— , TRUE,  means  that 
snapshot  output  is  enabled,  .FALSE,  means  snapshot  output  is  disabled. 

Simulation  Parameters  and  Results  Summary  Option  ( run_info ) 

The  simulation  parameters  and  results  summary  file  (run_id/ /' _ info.txt’)  is  optional,  but  is 
highly  recommended  as  it  contains  all  of  the  pertinent  inputs  and  results  of  the  simulation. 
run _  info— .TRUE,  means  that  summary  output  is  enabled,  .FALSE,  means  summary  output 
is  disabled. 

Exposure  Data  Matrix  Option  (exp^data) 

The  exposure  data  matrix  file  (run_id/ /'  _e dm5)  is  optional;  since,  it  is  needed  only  when 
examining  the  exposure  of  individual  elements  rather  than  being  concerned  with  the  overall 
exposure  of  the  antenna,  ex data—, TRUE,  means  that  the  exposure  data  matrix  output  is 
enabled,  .FALSE,  means  that  the  exposure  data  matrix  output  is  disabled. 

Temporal  Wave  Data  Option  (wave__temporal) 

The  temporal  wave  data  file  {run_id/ /' _ mat2’)  is  optional,  since  it  is  needed  only  when  an¬ 
alyzing  the  power  spectrum  of  the  modeled  seas  in  time.  It  contains  the  wave  elevation  at  a 
single  random  point  on  the  antenna  in  time.  wave_  temporal—, TRUE,  means  that  the  temporal 
wave  data  output  is  enabled,  .FALSE,  means  that  the  temporal  wave  data  output  is  disabled. 

Spatial  Wave  Data  Option  (wave _  spatial) 

The  spatial  wave  data  file  ( run_  id/ /’_mat3’)  is  optional,  since  it  is  needed  only  when  analyzing 
the  power  spectrum  of  modeled  seas  along  the  antenna.  It  contains  the  wave  elevation  at  every 
point  on  the  antenna  for  the  time  —  0.0  seconds.  wave_  spatial=. TRUE,  means  that  the  spatial 
wave  data  output  is  enabled,  .FALSE,  means  that  the  spatial  wave  data  output  is  disabled. 
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Significant  Wave  Height  ( sig_  wave_  ht ) 

This  variable  sets  the  sea’s  significant  wave  height,  a  parameter  used  to  describe  the  sea  severity 
in  both  the  Bretschneider  and  Ochi  spectral  models.  Appendix  B  correlates  sea  state  to 
significant  wave  height  for  a  fully  developed  sea. 

Development  ( development ) 

This  variable  allows  the  user  to  shift  the  sea  spectral  energy  modal  frequency  up,  for  developing 
seas,  or  down,  for  decaying  seas.  This  parameter  applies  only  to  the  Bretschneider  spectral 
model. 

Spectral  Model  (spectral _model) 

Four  spectral  model  subroutines  are  available:  the  two  parameter  Bretschneider  spectrum, 
the  Ochi  “most  probable”  six  parameter  spectrum,  a  square  wave  spectrum,  and  a  single  sine 
wave  spectrum.  Option’!’  means  Bretschneider,  ’2’  means  Ochi,  ’3’  means  square,  and  ’4’  means 
sinusoidal.  The  square  wave  spectrum  and  the  sine  wave  spectrum  were  included  for  debugging 
and  validation  of  the  results.  Parameters  of  the  square  wave  spectrum  cannot  be  changed  with 
the  input  file,  but  the  sine  wave  amplitude  and  frequency  are  computed  from  the  sig_  wave_  ht 
and  development  variables.  The  user  is  directed  to  review  the  source  code  for  the  sine  wave 
subroutine  to  determine  how  these  two  input  variables  control  the  sine  wave  parameters. 

Selection  of  the  spectral  model  has  a  measurable  impact  upon  the  exposure  statistics  com¬ 
puted  by  the  code.  There  is  no  litmus  test  that  the  user  can  apply  in  selecting  the  correct 
spectral  model,  but  the  different  models  are  better  suited  for  some  specific  conditions.  The 
assumptions  used  in  deriving  the  Bretschneider  spectrum  from  North  Atlantic  sea  data  were 
deep  water,  unlimited  fetch,  uni-directional  seas,  and  no  swell.  The  Ochi  spectrum  was  derived 
using  the  same  assumptions  except  that  it  allows  for  swell  from  a  distant  storm.  Based  on 
those  assumptions,  one  can  see  that  the  Bretschneider  spectrum  is  more  suited  for  open  ocean 
wind  driven  seas,  and  the  Ochi  spectrum  is  more  suited  for  open  ocean  wind  and  swell  driven 
seas. 
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Directional  Spreading  Function  (spreading _func) 

This  allows  the  user  to  select  the  directional  spreading  function  to  apply  to  the  power  spectra 
to  obtain  the  desired  wave  field.  Two  directional  spreading  functions  are  available:  cos2  9  and 
cos4 9.  Option  ’1’  means  cos2 9  and  option  ’2’  means  cos49. 

Number  of  Propagation  Directions  ( NN ) 

This  allows  the  user  to  select  the  number  of  discrete  unidirectional  waves  which  will  combine 
to  form  the  multidirectional  sea  surface.  A  larger  number  of  directions  will  result  in  more 
realistic  sea  spectra  but  at  the  expense  of  computation  time.  It  is  recommended  that  the  user 
experiment  with  the  number  of  propagation  directions  until  the  simulation  results  converge. 

Trapezoidal  Integration  Convergence  Factor  ( epsilon ) 

This  refers  to  the  convergence  factor  of  the  trapezoidal  integration  that  is  needed  to  calculate 
the  integral  of  the  directional  spreading  function.  A  typical  value  is  le-3. 

Antenna  Tow  Angle  (tow_  angle) 

Antenna  tow  angle  with  respect  to  the  primary  axis  of  propagation.  A  tow  angle  of  0  degrees 
represents  following  seas,  90  degrees  represents  beam  seas,  and  180  degrees  represents  head 
seas. 

Antenna  Diameter  (diam_  antenna) 

This  is  the  outer  diameter  of  the  hose-like  structure  encasing  the  antenna.  Element  performance 
(gain)  is  controlled  by  whether  this  outer  diameter  is  submerged;  the  diameter  of  the  inner 
conducting  slot  antenna  affects  the  rotational  invariance  of  the  element,  and  the  maximum 
theoretical  gain  of  a  fully  exposed  element. 

Water  Density  ( rho_water ) 

Water  density  in  which  the  antenna  operates. 
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Drogue  Size  (drog^area) 

The  drogue  frontal  area  (in  m2)  determines  the  tension  at  the  trailing  end  of  the  antenna, 
according  to 

T  =  \ PwaterU 2  drog  _area  CD_n  (A.8) 

A  typical  value  is  0.01  m2. 

Normal  Drag  Coefficient  ( CD__n ) 

This  represents  the  drag  coefficient  for  a  cylinder  in  cross  flow.  A  constant  value  of  1.0  has 
been  assumed. 

Tangential  Drag  Coefficient  (CZ?_f) 

This  is  the  skin  friction  associated  with  flow  along  the  long  axis  of  the  antenna.  A  relatively 
smooth  outer  surface  has  been  assumed  in  choosing  this  value,  which  has  been  assumed  constant. 
A  typical  value  is  0.0035. 

Antenna  Density  (rho_  antenna) 

The  average  antenna  density,  which  is  assumed  to  be  constant  over  the  entire  array.  This  is 
an  approximation  in  that  the  antenna  has  slightly  higher  density  at  the  actual  element,  and 
slightly  lower  density  between  elements. 

Antenna  Lateral  Stiffness  (flex_  stiffs  coeff) 

This  variable  represents  the  resistance  of  the  antenna  to  bending,  and  is  given  by 

flex_stiff_coeff  =  El  (A.9) 


where 

E  is  the  modulus  of  elasticity 
I  is  the  area  moment  of  inertia. 
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Figure  A-l:  Sample  Threshold  Antenna  Exposure  of  0.04 
Added  Mass  Coefficient  ( C_m ) 

The  added  mass  coefficient  for  a  cylinder  in  cross  flow  has  been  set  to  1.0. 

Maximum  Number  of  Samples  (max_num_  of _  samples) 

(r) 

This  refers  to  the  maximum  number  of  “snapshots”  taken  for  the  Matlabw  viewer,  or  the 
exposure  data  matrix.  max_  num_  of_  samples  is  used  to  dynamically  allocate  array  storage  for 
variables  which  are  sampled  periodically  in  time. 

Threshold  Antenna  Exposure  ( reqd_  ant_  exposure) 

As  mentioned  earlier,  the  antenna  gain  remains  steady  until  the  element  just  becomes  sub¬ 
merged.  In  order  to  capture  this  performance,  and  introduce  margin  for  capillary  waves  not 
modeled,  the  parameter  reqd_  ant_  exposure  has  been  defined.  The  reqd_  ant_  exposure  is  the 
height  above  the  water  which  must  be  maintained  by  the  antenna  to  be  considered  exposed. 
The  reqd_  ant_  exposure  is  given  as  a  fraction  of  the  antenna  diameter.  See  Figure  A-l  for  a 
diagram  representing  an  antenna  diameter  of  1.0  and  a  reqd_ant_  exposure  of  0.04. 
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Threshold  Number  of  Elements  ( num_  of_  elem_  thresh ) 

This  variable  allows  the  user  to  specify  the  lowest  number  of  exposed  antenna  elements  necessary 
to  achieve  the  desired  design  antenna  array  gain.  At  each  stat_  time,  the  number  of  exposed 
elements  is  compared  to  num_  of _elem_  thresh]  the  result  (a  “1”  when  the  threshold  is  met, 
“0”  otherwise)  is  found  in  the  summary  output  file,  run_  id//’_info.txt’,  along  with  any  other 
variables  tracked  over  the  course  of  the  simulation.  Currently,  the  variables  tracked  over  time 
are 


•  tow  point  tension  ( tow_tension(  )) 

•  number  of  elements  exposed  ( is_exposed(  ,  )) 

•  threshold  met  ( thresh_met :  “0”  or  “1”,  as  above) 

•  counterdetection  length  ( length_totally_exposed(  )). 
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A.2.2 


Sample  Input  File 


Input  File  for  proteus  Simulation  Program 

SI  units  throughout;  useful  conversions  follow 
1  inch  =  0.0254  m 
1  foot  =  0.3048  m 
1  knot  =  0.5144  m/s 


15M  ay_02 

run_id 

8  characters  long 

13 

nnf 

#  of  components.  N_sea=2**nnf 

60 

num_of_ele 

#  of  actual  slotted  antenna  elements 

0.5 

element_length 

physical  slotted  antenna  length  (m) 

8 

M 

#  of  mesh  points/element_length 

3.0 

U 

tow  speed  (use  POSITIVE  #)  (m/s) 

30.0 

duration 

total  simulated  run  time  (sec) 

.001 

t_  incr 

time  step  increment  used  for  antenna  calculations 

0.2 

pic_time 

snapshot  interval  for  Matlab  viewer  (sec) 

0.0 

t_first_pic 

time  of  first  snapshot  (sec) 

0.1 

stat  time 

statistics  interval  (sec) 

.FALSE. 

matlab_output 

TRUE:  A  file  is  created  that  contains  wave  and  antenna  data 

.TRUE. 

run_  info 

TRUE:  A  file  is  created  that  contains  a  summary  of  run  parameters  and  results 

.FALSE. 

exp_data 

TRUE:  A  file  is  created  that  contains  the  exposure  data  matrix 

.FALSE. 

wave  temporal 

TRUE:  A  file  is  created  that  contains  temporal  wave  elevation  at  a  single  point 

.FALSE. 

wave_spatial 

TRUE:  A  file  is  created  that  contains  spatial  wave  elevation  at  time=0 

3.0 

sig_  wave_  ht 

significant  wave  height  (m) 

1.0 

development 

modal  frequency  shift  factor  (used  only  for  spectral  model  1) 

1 

spectral_model 

l  =  Bretschneider;  2=Ochi;  3  =  square;  4  =  sine 

1 

spreading_  func 

directional  spreading  function:  1  —  cos~2;  2=cos'*4 

9 

NN 

Number  of  propagation  directions  (must  be  odd  number) 

le-3 

epsilon 

Convergence  factor  for  trapezoidal  integration 

180.0 

tow  _  angle 

Tow  angle  relative  to  primary  axis  of  propagation  (0-180  degrees)(0  following, 180  head) 

0.0762 

diam  _  antenna 

Antenna  outer  diameter  (meters) 

1025.9 

rho_  water 

density  of  seawater  (kg/m'S) 

0.01 

drog  area 

drogue  area  normal  to  fluid  flow  (m'2) 

1.0 

CD  _  n 

normal  flow  drag  coefficient  for  a  cylinder 

0.0035 

CD_t 

antenna/fluid  tangential  drag  coefficient 

408.0 

rho  antenna 

antenna  density  (kg/ra^S) 

100.0 

flex__stiff_coeff 

lateral  bending  stiffness  of  antenna  (Pa  m~4)  El 

1.0 

C_  m 

added  mass  coefficient  for  a  cylinder 

1000 

m  ax  _  num  _  of  _  samples 

upper  limit  on  #  of  samples  in  time  domain 

0.04 

reqd_  ant_exposure 

threshold  amount  of  vertical  exposure  necess  for  reception  (%  of  antenna  diameter) 

30 

num  of_  elem  thresh 

threshold  #  of  elements  required  to  achieve  necessary  gain 

A. 3  Output  Files 

A. 3.1  Snapshot  Output  File  (run_id/ /’ _ mat’) 
File  Description 

This  file  provides: 
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Figure  A-2:  Sample  Output  of  Snapshot  Output  File 

•  a  time  history  of  sea  elevation  along  the  antenna 

•  a  time  history  of  vertical,  longitudinal,  and  transverse  sea  velocities  along  the  antenna 

•  a  time  history  of  vertical  antenna  position 

(r) 

•  antenna  parameters  necessary  for  the  Matlabw  function  file. 

(R) 

This  file  is  intended  to  be  viewed  with  the  Matlab^  function  file,  proteus _  mat 

Sample  Output 

(r) 

A  sample  output  of  the  Matlati  function  file  proteus_  mat  is  shown  in  Figure  A-2 

A.3.2  Simulation  Parameters  and  Results  File  (run_id/ _ info.txt’) 

File  Description 

This  file  provides: 

•  a  complete  listing  of  antenna  parameters 


Run  ID:  24May-01 

Sea  Data  For  Time  0.00  sec 
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•  a  complete  listing  sea  spectra  parameters 


•  a  summary  of  exposure  and  antenna  performance  statistics. 


This  file  is  intended  to  be  viewed  with  any  standard  text  editor. 


Sample  Output 


Summary  of  Parameters  and  Results 


General  Information 


Run  ID  =  24May  01 

Parameter  and  Result  Summary  File  =  24May_01_info.txt 
Wave  and  Antenna  Data  File  =  24May _01_mat 
Start  Date/Time  =  05/12/2000  07:22:31 
End  Date/Time  =  05/12/2000  07:41:07 
Elapsed  Computational  Time  =  1115.56  seconds 
Simulated  Analysis  Time  =  30.00  seconds 

Antenna  Parameters 


Number  of  Antenna  Elements  =  60 

Length  of  Antenna  Element  =  0.50  meters 

Total  Antenna  Length  =  30.00  meters 

Antenna  Diameter  =  0.0762  meters 

Lateral  Bending  Stiffness  —  50.00  N-m~2 

Antenna  Density  =  408.00  kg/m  "3 

Drogue  Area  Normal  to  Fluid  Flow  =  0.0100  m~2 

Normal  Flow  Drag  Coefficient  for  a  Cylinder  =  1.000000 

Tangential  Drag  Coefficient  =  0.003500 

Coefficient  of  Added  Mass  for  a  Cylinder  =  1.000000 
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Number  of  Mesh  Points  Per  Element  -  8 


Tow  Speed  =  2.50  m/scc 
Tow  Angle  =  180.0  degrees 

Time  Increment  for  Antenna  Related  Calculations  =  0.00010  seconds 

Sea  Spectra  Parameters 


Spectral  Model  =  Bretschneider 

Directional  Spreading  Function  =  cos'- 2 

Development  Factor  =  1.00 

Significant  Wave  Height  =  3.00  meters 

Minimum  Cutoff  Frequency  =  0.35  radians/sec 

Maximum  Cutoff  Frequency  =  22.21  radians/sec 

Length  of  Seas  Modeled  =  512.0  meters 

Number  of  Spatial  Sea  Domain  Mesh  Points  =  2"  13 

Number  of  Propagation  Directions  —  3 

Seawater  Density  =  1025.90  kg/m  "3 

Time  Increment  for  Sea  Spectra  Related  Calculations  =  0.02125  seconds 
Number  of  Antenna  Time  Increments  Between  Sea  Spectra  Calculations  —  212 

Simulation  Results 


Maximum  Tow  Tension  =  116.4  N 

Minimum  Tow  Tension  =  24.3  N 
Stretched  Antenna  Length  =  30.2  meters 
Antenna  Strain  —  .005 

Statistics  Sampling  Interval  =  0.100  seconds 
Number  of  Samples  Observed  for  Statistics  =  301 
Average  Fraction  of  Antenna  Exposed  =  0.895 

Vertical  Antenna  Exposure  Necessary  to  Receive  =  4.0  percent  of  diameter 
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Number  of  Antenna  Elements  Required  to  be  Exposed  =  30/60 


Time  That  Number  of  Required  Elements  are  Exposed  =  100.0  percent 


Average  Detection  Length  =  1.95  meters 


A.3.3  Exposure  Data  Matrix  File  (run_zd//’_edm’) 
File  Description 

This  file  provides: 


•  the  raw  exposure  data  information  for  each  antenna  element  at  each  time  sampled 

•  the  spatial  position  of  all  elements. 


This  allows  one  to  analyze  statistics  across  time  or  space.  This  file  can  be  viewed  with  a 
(R) 

user  supplied  Mat  lab'  file  or  with  a  spreadsheet  application. 

Sample  Output 

time  spatial  position  /  exposure  status 


0.00000 

0.09900 

0.19800 

0.29700 

0.39600 


30.00000  29.50000  29.00000  28.50000  28.00000 

* 

• 

1 

1 

* 

1 

0 

• 

1 

0 

1 

1 

1 

• 

0 

1 

1 

1 

1 

• 

1 

1 

1 

0 

1 

. 

A.3.4  Temporal  Wave  Data  File  (run_  id/ /’_mat2’) 
File  Description 

This  file  provides: 
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Power  Spectra!  Density  Comparison  for  Data  Run  <17May-02> 


<V_  8 


Scaled  Output  Spectrum 
Input  Spectrum 


Area  of  Input  Spectrum  =  5.000000 

Area  of  Unsealed  Output  Spectrum  (Not  Shown)  =  0.056341 
L72k  Scaling  Factor  =  86.84 
Area  of  Scaled  Output  Spectrum  =  4.892339 
Number  of  data  points  per  segment  =  4096 
Number  of  non-o\«r1apping  data  segments  =  2, 

Total  Number  of  data  points  =  l 


0.6  0.8 
frequency  (Hz) 


Figure  A-3:  Sample  Output  of  Temporal  Wave  Data  Analysis 

•  a  time  history  of  the  wave  elevation  at  a  single  random  point  on  the  antenna 

(r) 

•  antenna  parameters  necessary  for  the  Matlahi  function  file. 

This  allows  one  to  analyze  the  power  spectrum  of  the  modeled  seas  in  time.  It  is  intended 

to  be  processed  with  the  FORTRAN  executable  file,  proteus _  mat2.  exe,  then  viewed  with  the 
(R) 

Matlab^  function  file,  proteus_  mat2. 

Sample  Output 

A  sample  output  of  the  Malian'  function  file  proteus_  mat2  is  shown  in  Figure  A-3. 

A. 3. 5  Spatial  Wave  Data  File  (run_  id//’_mat3’) 

File  Description 

This  file  provides: 

•  the  wave  elevation  at  every  point  on  the  antenna  at  time  =  0.0  seconds 

•  antenna  parameters  necessary  for  the  Mat  lab  7  function  file. 


96 


16 

14 

12  - 

10 


Power  Spectral  Density  Comparison  for  Data  Run  <17May-01> 


Scaled  Output  Spectrum 
Input  Spectrum 


Area  of  Input  Spectrum  =  5.000000 

Area  of  Unsealed  Output  Spectrum  (Not  Shown)  =  0.075695 

L/2k  Scaling  Factor  =  64.00 

Area  of  Scaled  Output  Spectrum  =  4.844453 

Number  of  data  points  per  segment  =  1024 

Number  of  non-overlapping  data  segments  =  8 

Total  Number  of  data  points  =  8192 

NOTE:  Output  spectrum  data  is  at  time  t  =  0 


spatial  frequency  (m  ) 


Figure  A-4:  Sample  Output  of  Spatial  Wave  Data  Analysis 


This  allows  one  to  analyze  the  power  spectrum  of  the  modeled  seas  in  space.  It  is  intended 
to  be  processed  with  the  FORTRAN  executable  file,  proteus _mat3.exe,  then  viewed  with  the 
Matlabw  function  file,  proteus_  mats. 


Sample  Output 


A  sample  output  of  the  Matlabw  function  file  proteus_  mats  is  shown  in  Figure  A-4. 
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Appendix  B 


Sea  State  Correlation  Data 
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Significant ' 

Modal 

Wind 

Sea 

Wave  Height 

Frequency 

Speed 

State 

(m) 

(rad/sec) 

(knots) 

0.1 

3.96 

4.2 

0.2 

2.80 

5.9 

0.3 

2.29 

7.3 

1 

0.4 

1.98 

8.4 

0.5 

1.77 

9.4 

0.6 

1.62 

10.3 

0.7 

1.50 

11.1 

0.8 

1.40 

11.9 

2 

0.9 

1.32 

12.6 

1.0 

1.25 

13.3 

1.1 

1.20 

13.9 

1.2 

1.14 

14.6 

1.3 

1.10 

15.2 

3 

1.4 

1.06 

15.7 

1.5 

1.02 

16.3 

1.6 

0.99 

16.8 

1.7 

0.96 

17.3 

1.8 

0.93 

17.8 

1.9 

0.91 

18.3 

2.0 

0.89 

18.8 

4 

2.2 

0.85 

19.7 

2.4 

0.81 

20.6 

2.6 

0.78 

21.4 

2.8 

0.75 

22.2 

3.0 

0.72 

23.0 

5 

3.2 

0.70 

23.8 

3.4 

0.68 

24.5 

3.6 

0.66 

25.2 

6 

3.8 

0.64 

25.9 

4.0 

0.63 

26.6 

4.5 

0.59 

28.2 

5.0 

0.56 

29.7 

5.5 

0.53 

31.2 

6.0 

0.51 

32.6 

6.5 

0.49 

33.9 

7 

7.0 

0.47 

35.2 

7.5 

0.46 

36.4 

8.0 

0.44 

37.6 

8.5 

0.43 

38.8 

9.0 

0.42 

39.9 

Table  B.l:  Sea  State  Correlation  Data 
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